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Abstract 

Coherent structures, such as solitary waves, appear in many physical problems, including fluid me- 
, chanics, optics, quantum physics, and plasma physics. A less studied setting is found in geophysics, where 

cl ■ 

highly viscous fluids couple to evolving material parameters to model partially molten rock, magma, in 
the Earth's interior. Solitary waves are also found here, but the equations lack useful mathematical 
structures such as an inverse scattering transform or even a variational formulation. 

A common question in all of these applications is whether or not these structures are stable to 

(N ■ 

^ ■ perturbation. We prove that the solitary waves in this Earth science setting are asymptotically stable and 

cn 

' extending beyond the physical parameter space. Furthermore, this extends existing results on well 



accomplish this without any pre-exisiting Lyapunov stability. This holds true for a family of equations, 
extending beyond the physical parameter space. Further 
posedness to data in a neighborhood of the solitary waves. 
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1 Introduction 



Coherent structures, such as solitary waves, appear in many physical problems, including fluid mechanics, 
optics, quantum physics, and plasma physics. A less studied setting is found in geophysics, where highly 
viscous fluids couple to evolving material parameters to model partially molten rock, magma, in the Earth's 
interior. Solitary waves are also found here, but the equations lack the useful structures such as an inverse 
scattering transform or even a variational formulation. 

A important question in all of these applications is whether or not these coherent structures are stable to 
perturbation. We prove that the solitary waves in this Earth science setting are asymptotically stable and 
accomplish this without any pre-exisiting Lyapunov stability. 

1.1 Magma— Porous Flow in a Viscously Deformable Media 

Models of magma in the Earth's interior couple Stokes flow of the viscous melt to the slow, creeping de- 
formation of the porous rock. These stress balance equations couple to transport equations for the volume 
fraction of melt, the porosity. Formulations may be found in [5, 19,34,35,40]. Nonlinearity appears in fluxes 
and through the material properties, the permeability and viscosity of the porous, deformable rock, which 
depend nonlinear ly on the porosity. Consequently, such models are known, from computations, to feature 
localization and generate coherent structures, see [1-3,15,34,35,40-43,50]. The physical assumptions and 
their implications will be discussed in a forthcoming review article, [37]. 

Under certain assumptions (small fluid fraction, absences of large-scale shear, no melting, etc.), such 
a system reduces to a single scalar equation for the porosity's evolution, [2,3,40,41]. The d-dimensional 
equation is 

M + ir) - V • [0"V (<^-™at0)] =0, X e t > (1.1) 

with the boundary conditions that 0(x, i) 1 as ]x] oo. V — {dx,dy,dz) for d = 3 and V ~ {dx,dz) 
for d — 2. The nonlinearity n comes from the relationship between the permeability, K, and the porosity of 
the rock, i^T cx 0". m relates to the bulk viscosity, C, to the porosity of the rock, C oc 0"™. In the physical 
regime, these exponents have values 2 < n < 3 and < m < 1, [12,13,19,30,34,35,44,45,51]. 

Equation (|l.ip appears elsewhere in Earth science as a model for convective Mantle plumes. Referred to 
as hot .spots at the surface, these plumes are localized regions of upwelling hot, buoyant material. Examples 
include the Hawaiian Island chain and Iceland. Modeled as the flow of a viscous fluid up a conduit embedded 
in a higher viscosity medium, an equation of the form (jl.ip was derived in [23]. There, the equation is one- 
dimensional {d — 1), the exponents are {n,m) — (2, 1), and the depdendent variable is the cross-sectional 
area of the pipe. 

Numerical simulations of (jl.ip . in one-, two-, and three-dimensions were performed in [2,3,34,35,50], 
where stable, radially symmetric, solitary traveling waves were observed. In [22], it was shown that in 
one-dimension, solitary waves, Uc{x — ct), in excess of the reference state, 0=1, exist if n > 1. In the 
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Figure 1: Solitary waves colliding and propagating in an experiment with honey. Figure 1 of [36]. 

context of conduit flow, discussed in the preceding paragraph, analog experiments using viscous syrups 
appear in [23,36,48]; robust solitary waves appeared as predicted, see Figure [TJ 

1.2 Stability of Solitary Waves 

We consider (jl.ip in one-dimension, 

dtcf> + d, (0") - 5, = 0, lim cj,{x, t) = 1 (1.2) 

\x\ — >oo 

where the z coordinate has been relabeled x. A cursory explanation for the solitary waves' stability may be 
found in [49]. Under a small amplitude scaling, (|1.2p is, to leading order, governed by the Korteweg-de Vries 
(KdV) equation. Since KdV solitons are stable, on a timescale for which KdV approximates (|1.2p its solitary 
waves should also be stable. 

Based on observations of numerical experiments, we expect a slightly perturbed solitary wave to evolve 
into another wave with similar amplitude and phase. It will be accompanied by some small amplitude 
dispersive waves and, perhaps, another solitary wave of smaller amplitude. The leading wave will outrun 
these other disturbances, cease interacting with them, and stabilize. 

Some intuition for this stability may be found in two properties. First, taller solitary waves travel with 
greater speed, c, than smaller ones. In the frame of the largest solitary wave, y — x ~ ct, the other waves 
travel leftwards. Second, in the frame of the leading solitary wave, small perturbations of the reference state, 
(j){x, t) — 1 + ■ip{x — ct, t) and ^ 1, evolve under the linear flow 

df-tp - cdyip + ndyip - dldtip + cdyip = 
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The dispersion relation and group velocity are 



uj^k) = — ck, Lu'{k) = n- 



Since both phase and group velocities are negative for all fc, small dispersive waves also travel leftwards. 
These two mechanisms are diagrammed in Figure [51 and we will exploit them to prove the main theorems. 

As the system is conservative, perturbations, such as a small solitary wave, will not vanish in a translation 
invariant norm. A suitable norm will register leftward motion as decay. We will use exponentially weighted 
norms, in the frame of the leading solitary wave. These norms are defined in Section [1.51 

Several paths to proving stability are available. One method is to seek constants of motion that can be 
combined into a metric centered at the solitary wave. Since the metric is time independent, if the perturbation 
is initially small, it will remain so. This elegant method relies on the calculus of variations and for equations 
such as KdV and Nonlinear Schrodinger (NLS) it may be used to prove orbital stability, [4,7,46,47]. A 
solitary wave, J7c, is said to be orbitally stable if for data sufficiently close to it 

inf ||it(t) - UA- + v)\\x < S for all t 

for some S > and an appropriate norm \\-\\x- Typically, the norm is equivalent to L'^(M) or H^{M.). 

However, in general, (jl.2p lacks a sufficient number of conservation laws for this approach. Indeed, in [3], 
the authors searched for an additional conservation law, in hopes of proving orbital stability, [8]. There 
are many such equations, including some of the Boussinesq systems, [28], and many of the "compacton" 
equations, [31], which also lack such structure, yet appear, in numerical experiments, to possess stable 
solitary waves. Note that when n + m = 0, (|1.2p is Hamiltonian, and we have investigated this, proving 
orbital stability in [39]. We wish to consider the general case, which includes the physically interesting cases 
(n,m) = (2,1) and (3,0). 

Another approach to stability is to linearize the problem Ut = N{u,Ux,Uxx: ■ ■ •) about a solitary wave 
and establish linear stability. Then, one seeks a way to perturbatively "boost" this to prove stability for the 
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nonlinear flow. This may rely on direct spectral analysis of the linearized evolution operator. We employ 
this method, following [20,27,28]. Through this, we prove that the solitary waves are asymptotically stable, 
our main result. By this we mean, that in an appropriate norm, 

\\u{t) - UcWy ^ ast^+cx) 

for data sufficiently close to Uc- 

We note that another method recently appeared in [f8]. Without linearizing, the authors employ a virial 
inequality to directly prove asymptotic stability of gKdV solitary waves. 

Our problem is an example of an equation for which one can prove asymptotic stability in the absence 
of orbital stability. Upon reflection, it is clear that the asymptotic stability of generalized KdV and BBM 
solitary waves could have been proven without using the orbital stability results. 

1.3 Main Results and Outline 

The main results are: 

Theorem 1.1 There exists d, > n such that for all cq G (n,c*], if (j)ca{x ~ 0o) is a solitary wave solution of 
(|1.2p . then there exist constants > 0, e* > and a > 0, such that for e < e*, if 

then 

(a) (|1.2p has a solution with data 4>o{x) ~ (pcoix + 0o) + ^^0(2;) for all time. 

(b) There exist Coo, 9 00, and k > such that 

- (f'c^i- - Coot + 9oo)\\h^ <K,e (1.3) 
lie"" Coot -^^oo,^)-0c^(•)]||ffl <K,ee-^' (1.4) 

|coo - Col + 1^00 - ^ol <^*e (1-5) 

Corollary 1.2 Let n + m = 0. If dcM[4)c] > 0, Af defined in p.8|) . then Theorem \l.l\ holds for all c > c*, 

except for a discrete set with no accumulation point. 

Remark 1.3 Theorem \1.1\ and Corollary may be extended for all n > 1 and c > c^, up to the acceptance 
of a time independent numerical computation. 

The feature of (|1.2p that allows us to prove nonlinear stability from the linear stability is a non-negative 
invariant, denoted Af [</)], and defined in (|2.8p . Taylor expanding J\f about a solitary wave, 

+ f] = AA[0,] + {SM[<j)c],v) + {S^M[<j)c]v, v)+0 (M^H^) 
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The first variation does not vanish and the second variation is not a positive definite quadratic form. However, 
in the frame of the solitary wave, the perturbation v is migrating to — oo. Therefore 

{6Af[(j)c] ,v) ^0 as t ^ +00 

The second variation may be decomposed as S'^Af[(f>c] = P + Q, P & positive quadratic form and Q — Q{x) 
a locahzed function. Then since the perturbation moves leftward 

{Pv, v) > {v, v) 

{Qv, v) ^ as t —>■ +00 

Asymptotically, 

\\v\\h^<KAU 

giving a Lyapunov type bound on the perturbation. 

However, more is needed to formalize this into a proof, notably a sense in which the perturbation recedes 
from the solitary wave. This is accomplished by analyzing the spectrum of the linearized evolution operator 
in a weighted space, in which the perturbation will decay. 

The plan of the proof is as follows 

(I) In Section [2] we review properties of (|1.2p and establish regularity properties of the solitary waves. 

(II) In Section [31 we prove that the linearized operator, Aa, has the property that there exists e > such 
that 

a (Aa) n > -e} = {0} 
and zero is an eigenvalue of algebraic multiplicity two. 

(III) In Sectional we prove 

\\w{t)\\Hr = \\e^-'wo\\H^<Ke-'"\\w,\\H^ 
for appropriate wo, K and b positive constants. 

(IV) In Section [51 we make several estimates, including a formalization of the Lyapunov bound. We also 
formulate equations for the speed and phase parameters of the solitary wave (c(t), 0{t)), coupling them 
to the infinite dimensional system for the perturbation. 

(V) In Section [6l we prove the main results, asymptotic stability and global existence of data near a solitary 
wave solution. 

Some remarks are made in Section [71 and additional details are located in the Appendices. 
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1.5 Notation 

Generic constants will typically be denoted by the capital letters K, M, and N, sometimes with tildes, 
overlines, or primes. Subscripts, such as Mj, may appear to indicate that M depends on 7. We avoid using 
C as a generic constant, as c appears throughout the paper as the speed parameter, and an operator C(A) 
appears in Section [S] 

Functions will typically live in spaces iJ'''(M) = W^'''(R), k a non-negative integer, the spaces of square 
integrable functions with square integrable (weak) derivatives up to order k. We will frequently omit writing 
M. The LP{M.) spaces will also appear, in particular and L°° . While we write for the norm of a 

function in , we only write ||/||p for the norm of a function in . 

We will be interested in functions in the exponentially weighted space, 

H^ = {u: e'^^w(x) G H^] 
for fc = 0, 1, 2, . . ., and a > with associated norm 

We also define the norm 

\\f\\H-nHl^\\f\\H- + \\f\\Hl 

The exponential weight will always be a positive number; we will often omit the assumption a > in 
statements. 

Frequently, we will have an operator T defined on a weighted space, -ff^, but wish to make computations 
in the unweighted space. To T we associate Ta = e^^Te""^, an operator on H''. For the differentiation 
operator, dx '-^ Da — dx — a. 

2 Preliminaries 

2.1 Properties of the Equation in a Weighted Space 

Much of the analysis involves studying (jl.2|) in an exponentially weighted space. We therefore state the 
following extension of the well-posedness results obtained in [38] for spaces: 
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Theorem 2.1 (Local Existence in Time & Continuous Dependence Upon Data) 
Given < a < 1, let (j)o{x) satsify, 

1100 - l||//inHi < i? < oo 
inf (f>o{x) > ao > 

X 

inf 00 (a;)" - a^M^y > Po > 



(2.1) 
(2.2) 

(2.3) 



Then there exists Tiocai > and (p{x,t) — 1 G C^([0, Tiocai), H H^), a solution of (|1.2p with data 4>q, 
satisfying 



m-,t)~l\\HinHl <2i? 
inf (/)(x,t) > \ao 

M(P{x,t)"' -a^(pix,tr > IPo 



for t < Tiocal- 

Moreover, there is a maximal time of existence Tmax, such that if T-^^ax < oo, then 



(2.4) 
(2.5) 
(2.6) 



^Jim ||(/)(-,i) - lll^inifi + 



1 




1 




0(-,i) 


+ 

OO 


0(-,t)™ - a^(j){-,tY 


OC 



oo 



(2.7) 



Remark 2.2 When a = Q, (|2.3p is unnecessary; this case was treated in [38]. The importance of this 
condition for a > will be discussed in Section[ 



Theorem 2.3 Given < a < 1, ^^J) - 1 G C^{[0, T]; n H^), j = 1, 2 be two solutions of such that 

U'-^\;t)-l\\HinHl<B^<^ 

inf 0(J) {x, t) > ao > 

X 

inf(0(a;,i)(^))" - a\(j){x , t)'^^^^ > /3o > 
There exists a constant K = K{R^ ao^ [Sq^ a)^ such that 

||0W(.,i)-0('^(-,i)||ffinffl <e^*ll0^'^-0^"llHinff^ V ^ < 
Additionally, (II. 2|) possesses the conservation law 



/(ir'"02+01og(0)-0+l)da: 
/(ir'"02+0~l~log(0))dx 



if n + m = 1, 

if 71 + m = 2, 
i+(ra+m-2)(0-i) ^jj^ other n and 



(2.; 



{n+m — 1) (n+m — 2) 

Af is well defined for bounded from below away from zero and ||0 — < oo. It is also locally convex 



about 6=1. See Section 3 of 



for details. 
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2.2 Solitary Waves and their Analytic Properties 



Let us review the properties of the sohtary waves associated with p.2p . In particular, we identify their decay 
and regularity properties, and introduce the KdV scaling for later use. 

Substituting the traveling wave ansatz, (j)c{x,t) — (pdx — ct), into p.2p with boundary conditions 



y — >±oo 



we have, after one integration, 



lim 0c(y) = 1, lim 9^0c(2/) = for j = 1,2,. 



y^±oo 



Letting Wf, = 1 — (/)cj satisfies 



+ (Ue + 1)" - 1 + C [U, + 1)" dy ((Ue + 1)"" dyU,) = 



Equation (|2.10p may also be integrated up to a first order equation, 

ic"" (9.0c)'-i^i(</>c;c) = O 

after applying the boundary condition 0^ ~* 1 at ±oo. Fi depends on the particular exponents; 

+ (1 _ 1) £11^ _ 1 log(a;) ifm = l. 



Fi(x;c) = < 



^-l-(l-^)log(^)-^^ 
log(^) + (l-^)(i-l)-^^ 



'1 - 



1\ 



if n + TO = 1, 
if n + TO = 2, 
- otherwise. 



2— Ti— m V cl 1— n— m c 1 — m 

Using (j2.12p . an equivalent second order, self-adjoint, equation for the solitary waves is 

F2{x;c) = x™-" [-~{x - 1) +c^i (x" - 1) - 2mx"-'"'^^Fi{x;c)] 
Let us introduce the KdV scaling. Define 
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c — n 



Applying the scalings, 
(|2.1ip becomes 



^ = j{x-ct), Uc{y) = 7t/(?(2/);7) 

71 — 1 



1, 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 
(2.15) 



(2.16) 



(2.17) 



(2.18) 



Remark 2.4 The parameter 7, (|2.16p . will be used throughout the paper. Because it uniquely maps c € 
(n,oo) onto (0, 1), we will use c and 7 interchangeably. 
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We summarize what is known about (|2.1ip and (|2.18p in the following two results: 

Theorem 2.5 For any c > n > I, (|2.11|) has a unique positive, even solution Uc, going to zero at ±00. In 
the KdV scaling, (|2.18p . U is real analytic in the arguments (^,7) G K x [0, 1). When 7 = 



.2 /I 



f/(f;0) = [/40 =3sech^ (^-f 

Furthermore, for 7 in any compact subset 0/ [0, 1) 

9|f/(e; 7)e±« (sign(O)'' -> Kj{j) as ^ ^ ±00 for j = 0, 1, 2 

Corollary 2.6 Given a compact interval [0,7o] C [0, 1), there exists a constant K such that for all^ G [0:7o]; 

|9|C/(C;7)| < i^e-l^l /or j = 0,1,2, -00 < C < 00 (2.19) 

\dluc{y)\ < K-^ — -e-'>\y\ for j = 0, 1, 2, -cx3 < y < cx) (2.20) 

Proof: From [22], solitary waves exist provided c > n > 1 and m e R. Writing the problem as a two- 
dimensional dynamical system, we may apply the Stable Manifold Theorem about the hyperbolic point (0, 0) 
to deduce the exponential decay, as in Theorem 2.1 and Corollary 2.2 of [28]. 



Remark 2.7 When the parameter 7 is small, (pc is in the regime of small amplitude, long waves, where 
KdV appears as the leading order equation in a perturbation expansion of ()1.2|) . as in [49]. 

Corollary 2.8 For each c > n > 1, the solitary wave solution, (pc — 1, H^s in H'^{K). Furthermore, there 
exists (To > such that the solitary wave (j)c may be analytically continued off the real axis into the strip 
{z : |5z| < o-q}- 

Proof: This is a consequence of Corollarv l2.6l and Corollary 4.1.6 of [9], see Appendix lA.il 

■ 

Corollary 2.9 Given a solitary wave (j)c, c> n, assume < a < 7. Then (/)c — 1 G . 
Corollary 2.10 Let n> 1. 

(a) The mapping c^~^ (j)c ~ ^ is ((n, 00); i?'^) . In fact the mapping is analytic. 

(b) This mapping is analytic, and, for fixed x, a-^ 4>c{x) is analytic function of c. 
(b) Given a < i, the mapping is also ((n/ (1 — 4a^), cx)) ; n i?^) . 

Proof: All parts are proved using the implicit function theorem, applied to the functional 

T[c,u] = dlu + F2il + u;c) 

See Appendix IA.2I for details and [6] for a statement and proof of the implicit function theorem for analytic 
mappings. 
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2.3 Remarks and Assumptions on the Exponential Weight 

We see in Theorem 12.11 and Corollary 12.101 that the particular exponential weight restricts what data and 
which solitary waves will be permissible. For the solitary wave result, this restriction comes from the decay 
rate associated with the speed; see Corollarv 12.61 

In the case of the existence theorem, (|1.2I) may be written as 

^ (2.21) 

The operator is 

= (l)ixr ~ (2.22) 

This is a bounded operator on — s- provided (f> is continuous and bounded from below away from zero. 
However, the exponential weight introduces a second constraint. Consider solving -ff^u — f, f G for 
u £ H^. Letting g — e'^^f and v = e°^u, this is equivalent to solving 

[c/)™ - Da irDa)] v^g, veH^ 

Multiplying by v and integrating by parts, 

^™ _ a2^«) „2 _^ ^Q^^^2 dx^ j gvdx 

A unique solution exists, provided a and satisfy inf^; (pi^x)"^ — a^(/)(a;)" > 0; this is condition (|2.3p . 

We invert these restrictions; given a solitary wave of speed c, we will assume that a is sufficiently small 
so that these, and other, properties hold. There are three restrictions in what follows. 

Let 

ai = i7(c) (2.23) 
0c ~ 1 will then be in H^^ and , as will all solitary waves of nearby speed. Let 

02 = iinf0c(x)("-")/2 (2.24) 

Then 0^ will satisfy (|2.3p . with 

inf 0e(a;)"' - a^4>c{xT > ? inf 0c(x)" > 

Hence, the solitary waves will live in a set on which the existence theorem applies. Moreover, for all data 0o 
sufficiently close to 0c in the norm, an analogous lower bound will exist. 

Remark 2.11 inf^; 0c(a^)'™ "''^^ is related to a physical length scale known as the compaction length, [19J. 
This length, ^comp.; is given by 



It measures the distances over which there will be geometrical rearrangement of the material, appearing 
macroscopically as changes in (j), in response to viscous stresses. 
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The stipulation a < inf^ (t){xY"^ "^/^ may be interpreted as requiring the length scale associated with the 
exponential weight, a^^ , to never be smaller than this intrinsic, spatially varying, length. 

Finally, there is a constraints related to the essential spectrum of a linear operator, discussed in Section 
H Let 

^3 = ^ / 2c Z (2.25) 

]] n + 2c+ {n + 8c) V c 

This will ensure that for any a < 03, the essential spectrum is located in a specific part of the complex plane. 
Let 

a^(c) = min|ai(c),a2(c),ia3(c)| (2.26) 
Then for any a < a^,{c), all of these properties will be satisfied for ip sufficiently close in to 0^- 

2.4 Ansatz and Linearization 

Given a perturbed solitary wave solution, 0, of (|1.2p . assume that there exists decomposition of into a 
(time dependent) solitary wave of some speed c{t) and phase 9{t) and a perturbation, v; this decomposition's 
existence will be proved in Section FS.ll Transforming our coordinate system into the frame of this modulating 
solitary wave, 

y{x,t)^x- [ c{s)ds + 9{t) (2.27) 
"'0 

(l){x,t)^cj),(^f^{y{x,t))+v{y{x,t),t)^My,t) + v{y,t) (2.28) 
The perturbation, v, is governed by 

dtv = A,v - edyv - cd.cj), - edy(t>, + Ti[v; (2.29) 

where 

A,v = ^"^H-^^dyL.v (2.30) 
Lcv ^ -c^'^dl (C™«) + [c - n<j)-^ + cn{<j)-^ ~ 1)] v (2.31) 

and (|B.6[) gives the defintion of J-i[v;(j)c]j composed of terms nonlinear in v. We make two remarks about 
(|2.29p . First, the linear operator Ac, is time dependent; c — c{t), and we would prefer to work with a time 
independent linear operator. Second, the appearance of the term 9dyV will prove problematic to studying 
the equation in H^. 

The first problem is addressed by adding and subtracting Ac,, , and considering the difference Ac — Ac„ 
as another perturbation of the linear flow. To remove the dyV term, we introduce a renormalized time, 



c{s)ds - e{t) 







(2.32) 



The asymptotic stability proof of BBM required a similar transformation, [20] . In addition, the problem will 
be considered in the weighted space i?^, with w{y,t) = e"^w(2/, i). All together, we have: 
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Proposition 2.12 (Perturbation Equation) 

The perturbation to a solitary wave of speed Co associated with the ansatz in (j2.28p . v and its weighted 
perturbation w = e'^^v, evolves according to the equations in t and T-time respectively. 



dtv = AcV - edyv - cdcctyc - (^dy^c + ^\ b; 0c] 

dtw = Ac^aW -6{dy~a)w- e"^ {cdc(t)c + Sdycjjc^ + Gi[w, v; 0c] 

drv = Aco ^ (cdc(t>c + Ody(j)c] + S'[co, c, e]v + -^^i [v\ 0c] 

r ~ H ^ / r — H 



drW = Ac„nW 



.'^y (cdc^c + Ody(t>,) + Sa[co, c, 9]w + -^Gilw, V- 0c] 



= Ac„nW + Q 



w,ti;co,c, I 



(2.33) 
(2.34) 

(2.35) 
(2.36) 



The operator S and the terms Ti and Qi are given explicitly in AvvendixWl 

Proof: Qi is obtained from by substituting e~°^ui for one of the v^s; e'^^Ti [u; 0c] = Gi [e"^?;, v\ 0c]. The 
details appear in Appendix iBl 



N.B. From here on, we assume cq to be fixed and will suppress its appearance in the linear operators Ac^ 
and Aco,a- 



3 Spectral Properties of Linearization about a Solitary Wave 

In this section we analyze the spectrum of A and Aa- We will use c in place of cq and x in place of y as the 
independent variable. 

AY = Ay (3.1) 
A^{l~d., [red. (07"-)] (3.2) 
Lc = -c^dl (0-'"-) + [c - n0-i + cn (0-^ - l)] (3.3) 
Aa = e^'^Ae-"'' (3.4) 

Our goal is to prove that the linearized problem is asymptotically stable; ||e"^°*ti;ollHi — > as i ^ +oo. We 
will actually show something much stronger, that this convergence to zero happens exponentially fast. Our 
strategy is that of [20,28]. We will 

1. Identify the essential spectrum of Aa by showing it to be a compact perturbation of a constant coefficient 
operator. 

2. Rule out point spectrum (eigenvalues of finite multiplicity) of Aa for |A| sufficiently large via an operator 
estimate. 
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Figure 3: The spectrum of the operator Aa- The only eigenvalue with 5RA > — e is A = 0. 

3. Use the Evans function, an infinite dimensional analog of the characteristic polynomial, to rule out 
nonzero point spectra of Aa in the set of "small" A, which will be compact. 

4. Show decay in time of the Co-semigroup e"*"* associated with Aa- 

The spectral analysis is handled this section, and the semigroup theory in the following section. 
The principle result of this section is: 

Theorem 3.1 (Spectrum of Linearized Operator) 

(a) Let a € (0,0^,(7)]. The essential spectrum of Aa denoted by aessiAa) is a curve lying in the open left 

half-plane, with rightmost point —oj, 

— bj = max{3?z|z £ a^ss{Aa)\ < (3.5) 

(b) There exists 7^, G (0,1) such that for each 7 G (0,7*] and a G (0,0^,(7)], there exists £(7, a) > such 

that the only eigenvalue of Aa with SRA > — e is X — and this is an eigenvalue of algebraic multiplicity 
two. 

(c) In the Hamiltonian case, n + m = 0, part (ii) may extended for 7 G (7*, 1) to all but a discrete set with 

no accumulation point. 

The spectrum is pictured in Figure\^ 
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3.1 Essential Specturm 

We make use of the definition of the Essential Spectrum of an operator from [32,33], that for a closed, 
densely defined operator A on a Banach space X , 

aess(A)- fl a{A + C) (3.6) 

where JC{X) denotes the set of compact operators on X . Other definitions are possible and well known; see 
Chapter IX of [11] for a discussion of how these definitions relate to one another. 

(j{A) \ acss{A) then consists of point spectra. This is so because, by Theorem 7.27 of [33], if A is not in 
the essential spectrum, then XI — A is Fredholm with index zero. Hence, it has closed range and a finite 
kernel. Therefore, it must be an eigenvalue of finite multiplicity. 

To prove Theorem l3.1l Dart (i), we express Aa as a compact perturbation of a constant coefficient operator, 
A^, obtained by setting 4>c{x) equal to its asymptotic state, 1. 

A^ = {I-Dl)-'D,{-cDl + c-n) (3.7) 

The difference between Aa and A^ , given explicitly in (|B.16[) . may be shown to be an -compact 
operator. Hence, Aa is a compact perturbation of and 

{Aa) {AT) 

Upon examination of its Fourier symbol of , the essential spectrum of Aa is 



[il- a){-c{it- a)^ + c-n) 

-ess{Aa)={ 1_(,^_,)2 ' 



(3.8) 



and 



OjTl 

-uj = max{5Rz|z e aess{Aa)\ = -ac+ < 

1 — a"^ 



and aess{Aa) lies in the open left half plane if < a < 7. crcss(Aa) is pictured in Figured 

In addition, for < a < a^, < 03, 03 defined by (j2.25p . then the spectrum moves rightward as a ^ 0. This 
is because 03 is the value for which —lo is leftmost in C, maximizing the rate of decay of e"** as y ^ —00. 

3.2 Large Eigenvalues 

As in [20], we will study the eigenvalues of Aa by considering separately a large |A| regime and a small |A| 
regime. 

Rewriting the linear operator Aa as 

Aa = cC^a (C™-) - (C'O + ^^^^^"1 (C'^.^c') + Cn4>^ H^l^D a [(C' " l) '] (3-9) 

we note that A is an I? eigenvalue of Aa if and only if it is also an 1? eigenvalue of Aa — 0^7™ Aqi/)™, given 

by 

Aa = cDa - nH-^l^Da (C^'O + '^^ii-^la (C^'^.^/.,-) + CuH-la^ a [C (C' " l) '] (3.10) 
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We will rule out eigenvalues of Aa, thus ruling them out for Aa- This is equivalent to studying Aa in a space 
weighted by (j)~"^{x), a strictly positive, smooth, and bounded function. Aa is also a compact perturbation 
of A'^ = A^ ; they share the same essential spectrum. 
Let the operator C(A) satisfy 

Proposition 3.2 (a) The operator C(A) is compact for X not in <Tess- In particular, C(A) is compact for 
all X with 3fiA > —uj. 

(b) For any A G C \ acss{Aa), we have that X is an eigenvalue of Aa if and only if 1 E cr{C{X)) . 

(c) Let X G C\<Tcss(^a)- A sufficient condition for X not to he an eigenvalue of Aa is that ||C(A)|| < 1, with 

norm either L? or depending on which space is under consideration. 

Proof: Using the equivalence of eigenvalues of Aa and Aa, parts (b) and (c) will follow once (a) is 
established, see [20]. The Fourier symbol of {XI — A'^)^^ is 

A(l - {i£ - a)2) - {i£ - a){-cii£ -a)^+c-n) 

This operator is bounded for A not in the essential spectrum. The difference, given explicitly in (|B.17[) .is a 
sum of Hilbert-Schmidt compact operators composed with bounded operators on L^, hence C(A) is 

compact on this space. 

For C(A) to be a compact operator on H^, it will be sufficient to prove 

(/ - diy/'C{X){I - dlr'/^ = {XI AT)-\I ~ SilY''\Aa - A^){I - 

is compact on . {XI — A^)^^ is still bounded, and by commuting operators, it may be proven that 

(/ - 92)i/2(i^ _ Af){I - a2)-i/2 is compact. 

■ 

We now rule out eigenvalues outside a rectangle that scales with 7^. 
Proposition 3.3 Let 5 £ (0, 1) be fixed. 

(a) For any > n, there exists an M > such that for n < c < d,, if 

5RA > -^acj^ 

and either 

|!3A| > cMj^ or 5RA > cM-f^ 

then |lC(A)||L2^i2 <1-S. 

(b) This result also holds for \\C'{X)\\hi^h^ ■ 
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Proof; If AI > 1, then by the assumptions on SRA and 3A, A is not in (Tess(^a), hence we may apply 
Proposition 13.21 (a) to conchide C(A) is a compact operator. If the norm of C(A) is less than one, part (c) 
of that proposition will imply it is not an eigenvalue; we seek an M > 1 for which ||C(A)|| can be made 
sufficiently small. 

For all (7, a) in {(7, a) | 7 e [0, 70], a < 0^(7)}, there exist Ki and K2 such that 

\\C{X)h2_,L2 < \\{XI ~ A^y' H^}Da\\L^-.L^ (cK^j^) + \\{XI - A^y' H-lU^^L-^ {cK^i^) 

This comes from expanding the difference Aa—A"^ and commutating operators, see (jB.17p . (A/ — A'^) ^ H^^Da 
and {XI - A^y H~l are constant coefficient operators and we will treat them in Fourier space. 
Thus, if we can prove that for A satisfying the the hypotheses 



sup 



sup 



{i£ - a) cKi-f^ 



af - {ti - a) -c {t£ - af + c^i"^ 



< 



A 



l-it£- af - {ti - a) -c {t£ - af + cj^ 



< 



1-S 
\-6 



we will be done. 

Introducing the scalings A = cA7'^, I — 7^, and a = jd, this is equivalent to identifying M > 1 such that 
when 5RA > and either 5RA > M or |3A| > M then both 



sup 



sup 





1-6 


A 


1 - 72 (i^ - §f 


- {< - I?) 
K2 




< 2 
1-6 


A 




- « - -D) 




< 2 



(3.11) 
(3.12) 



are satisfied. 

Squaring both sides, p. lip and (j3.12p may be rewritten as two polynomial inequalities. Pi (5A, 5ftA, 1^) > 
and P2(3A, 5RA, ^) > 0, respectively. We will show that for appropriately chosen A, the inequalities hold for 
all ^. Pi and P2 are treated similarly. We study P2: 

P2 (5A, 5RA, = a(3A)2 + p^A + ryi {^Af + m^A + m 

a = + 2ei\i + 7'^') + (1 - i^e^f 

(3 = -2^7^ - 2^(1 + 7^(1 + 29'^)) - 2^(1 + (-3 + 7^)02 ^ ^2^4^ 
ryi = a 

772 = 26'^V + 26*^^(3 + 7^(-l + 26|2)) + 261(1 - e'^){l - -fO^) 
j]3=f + 9^{1 - e^f + ^"^(2 + W^) + C^(l + 36**) - (3^ 
13 = 2X2/(1 - 6) 



18 



Using similar analysis as for Pi in [20], we first consider P2 as quadratic in 5A. Examining its discriminant, 

discriminant = 7** 

- 36^^0^ + 0{f) 

If > 5RA > —0/2, then there exists Mq > such that for all QA, 7 G [0,70] and |^| > Mq the discriminant 
is negative and P2 > 0. Furthermore, since the coefficient a > 0, there exists i?i > such that P2 > when 
7 e [0,70], |5A| > > 3fiA > -61/2, and \C\ < Mo 

For7 e [0,70] and^ e (0,1), both r/i and 772 are positive; if Fzl^A, 5RA, > then Pal^A, 5RA + X, > 
for any K > 0. Therefore, P2 > for all C if 7 e [0,7o], |^^A| > i?i and > -0/2. Also, P2 > for 
Id > A/o and aU 3A if 3?A > -61/2. 

We must still treat the case of |5A| < Ri and < Mq simultaneously. Consider P2 as quadratic in 5RA. 
771 and 772 are positive for all ^ and 773 is bounded from below. Therefore there is some i?2 > such that 
P2 > for aU C if 7 e [0,70], |5A| < and 5RA > i?2. Thus, for any 6 e (0, 1), 70 e [0,7o], there exist 
i?i > and i?2 > such that P2 > for all ^ if 

7 G [0, 7o] 
SRA > -61/2 
|3A| < Ri or 5RA > R2 

For C(A) : i/^ H^, the proof is similar, with constants Ki and K2 in place of ii'i and K2. 



3.3 Small Eigenvalues— The Evans Function 

In this section we rule out eigenvalues of Aa in the set |A| < Mj^. This is done using the Evans function, 
an analytic function that vanishes at eigenvalues of Aa- The Evans function, D = D{X;j), is constructed 
for the eigenvalue problem using particular solutions of an associated dynamical system 

y = P(a:,A,7)y (3.13) 

y = 0(e^i"^) asx^+00 (3.14) 

and the adjoint system, 

i^-zB{x,X,j) (3.15) 

z = 0(e^i^) as a; ^-00 (3.16) 

fii will be the eigenvalue of smallest real part of the B°^\ the limit as x — > ±00 of B. When certain conditions, 
described in Theorem 13.41 are met, the Evans function exists and may be explicitly defined as 

I?(A;7)=z(a;;A,7)-y(x;A,7) (3.17) 
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The idea is to measure the angle between the subspace of solutions decaying at +00 with the subspace 
decaying at — cxd; hence the appearance of the dot product. The Evans has an equivalent formualation in 
terms of the determinant of the fundamental solution of (j3.13p . For a more complete discussion of the Evans 
function, see [26]. 

Theorem 3.4 [26, 28] 

Let VI he a simply connected subset ofC'^ . Suppose that the system (j3.13p satisfies the following hypotheses: 

(i) B : R X n ^ £riy.n continuous in {x, A, 7) and analytic in (A, 7) for fixed x. 

(ii) _B°°(A,7) = lim^^ioo S(a;, A, 7) exists for all (A, 7) G fi. The limit is attained uniformly on compact 

subsets of Q,. 

(iii) The integral 

/"OO 

||i?(x,A,7)-i?°°(A,7)||rfx 



converges for all (A, 7) € and the convergence is uniform on compact subsets ofVL. 

(iv) For every (A, 7) G il, the matrix (X^ j) has a unique eigenvalue of smallest real part, which is simple, 
denoted /ii. 

Then D{X; 7) is well defined and analytic on fl, such that D{\; ^) — if and only if \3.13\ has a solution y{x) 
satisfying (|3.14p and 

y(x)=o(e^i^) as x-> -00 (3.18) 

3.3.1 The KdV Evans Function 

In the case of the KdV equation, the eigenvalue problem may be scaled to 

d^LKdvY = (y-dlY + Y - 3scch' Qx^ = KY (3.19) 

Because the speed parameter has been scaled out, there is only one eigenvalue parameter, A. 
Making the identification 



y satisfies the the dynamical system 



y=(r, d^Y, LKdvYf (3.20) 



y-SKdv(a^,A)y (3.21) 
/ 1 \ 



1 -3sech(ia;)^ -1 



A y 

A complete description of the associated Evans may be found in [27]. We summarize: 
Theorem 3.5 (The KdV Evans Function ) 
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(3.22) 



(a) The Evans function Z?Kdv(A) associated with (I3.19P is given by 



i?Kdv(A) = 



Mi(A) + l 
Mi(A)-l 



where /ii(A) denotes the root o/ /i^ — /i + A = o/ minimal real part. 

(b) T/ie domain o/-DKdv(A) is t/ie sto complex plane 

AKdv = C \ ^-oo, 

(c) T/ie essential spectrum of AxdV ■ ^ curve contained entirely in the domain {A : 9?A < — e} 

/or some e > 0. Furthermore, if A^^y{a) denotes the component o/ C \ cress(^Kdv) that contains the 
right half-plane, then _DKdv(A) has no zeros in A^^y{a) except for a zero of multiplicity two at A — 0. 



3.3.2 The Evans Function applied 

The eigenvalue problem AY = XY, 

{I - [r.d., (</>-"•)] Y' 5- {-cC^I (C"-) + [c - + cn (C' - 1)] } y = AF 
is equivalent to 



d^L.Y - A [/ - a, (0c""-))] Y = 



Defining 



y = 



'Y, L.y + Aca^c^^^))^ 



y solves the dynamical system 



y = A,c)y 
I 



B{x, A, c) 



V 



-1-1)] A/c -c-Vc"" 




\ 



The matrix B may be decomposed as S = B°°(X, c) + R{x, A, c). 



i?-(A,c) = 



R{x,X,j) 



1 

7^ A/c -c-i 
A 









A (<^™ - 1) 



(3.23) 
(3.24) 

(3.25) 
(3.26) 



(3.27) 



(3.28) 
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Also note that for the corresponding adjoint eigenvalue problem 



L*d,W - A [/ - i€d.-)] W = 



under the identifications 



z solves 



z = —zB{x, A, 7) 
Theorem 3.6 (Properties of the Evans Function) 
(a) The Evans function is defined and analytic on the set C C^, 



(3.29) 



(3.30) 



(3.31) 



= {(-^, 7) I 7 e (0, 1) and A £ n^} 



with 



where 



n^^{X\m> -Ao} \ (-Ao, -Aeut(7)] 



(3.32) 



(3.33) 



-\/ 8c2 + 20cn - ^2 - ScVn^ + 8cn - n\/n'^ + 8c 



zn-f^ H pn7^ + 0(7'^) 



(3.34) 



3V3 9\/3 
and 

(b) Given (A, 7) G il, and a < a*(7). If X is to the right of acss{Aa) , then the following are equivalent 

• ^(A;7) = 

• X is an eigenvalue of Aa 

(c) For such zeros of D, the algebraic multiplicity of X as an eigenvalue of Aa is equal to the order of X as 

a zero of D{X; 7). 

(d) D(0;7) = d\D{0;^) = 0, hence it is an eigenvalue of algebraic multiplicity at least two. 



Remark 3.7 With this construction, for^X < 0, A G fly, the characteristic polynomial, p.35|) not only has 
a unique root of minimal real part, all roots have distinct real part. This is stronger than is needed. 

Remark 3.8 Acut is labeled as such because there is a branch cut in the Evans function there. 

Proof : Part (a) requires the verification of the hypotheses of Theorem 13.41 for this system. Applying the 
properties of the (j)c, Corollaries 12.81 and 12. 101 and examination of p.26p . B is clearly continuous in its three 
arguments for A G C and c > n. In addition, for fixed x, it will be analytic in (A, c), or, equivalently, (A, 7). 
Thus property (i) holds. 
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By Corollarv l2.61 the limiting matrix exists and B — B°° is in . This will be uniform on compact 
subsets of C X [0, 1), establishing properties (ii) and (iii). 

Lastly, we must verify the existence of ^i, the unique eigenvalue of minimal real part. We divide this 
into two parts, 5ftA > and 5RA < 0. The characteristic polynomial, Vifi), of B°° is 



cP{n) = {X- c^i)il- H^) (3.35) 



Following the analysis in Section 2(c) of [26] for a similar polynomial in the case of gBBM, one confirms that 
property (iv) holds for 5RA > and all 7, hence the Evans function exists in {5RA > 0} x [0, 1). 

Using the analysis in [20] for Theorem 2.7 of the polynomial, one can conclude the existence of some 
Ai > and identify ^^(7), such that for all 7 G (0, 1), a unique root of minimal real part exists for A in the 
set 

{A: 3fiA<-Ai}\(-Ai,-f2(7)] 

Alternatively, we give a more precise analysis of the (|3.35p in Appendix [C] that yields values of Aq and f2 
given in the proposition. This concludes the proof of part (a) . 

To prove part (b), we need a lemma regarding the location of f7css(^a)- 

Lemma 3.9 Let 7 e (0,1) and a < a*(7). Let ~ r2+(7,a) denote the component of C \ cress(^a) 
containing the origin. Then for A e fi+ H fly, the roots of the characteristic polynomial satisfy the relation 

^fii <-a< ^Hj^i (3.36) 

Proof: By inspection, if A <E CTcss(^a), there is a root of ()3.35|) with 5ft/ij = —a. Conversely, if there is a 
root with real part —a, then A is in the essential spectrum. Hence the characteristic polynomial has a root 
with real part —a, if and only if A e (Tcss(^a) 

As noted in [26], Section 2 (c), for large |A|, the roots of the characteristic polynomial (|3.35p are 

-l + 0(|A|-i), l + 0(|A|-i), A/c + 0(|A|-i) 

So for large A in the right half plane, (|3.36[) holds because a < 1. Now suppose for some A G il+ D fl~f the 
inequality were false. Because the ^fij depend continuously on A, equality would have to hold for some A, 
but then it must be that A G (Tcss(^a), which we have assumed is not the case. 



We now prove part (b) of Theorem 13.61 If D{X;j) — 0, then there is a solution to the ODE, y = By 
such that 

y(a;) 0{e'^^'^) as a; ^ +00 and y(a;) = 0(6^^"^) as a; — > -00 
Hence, Part 1(d) of Proposition 1.6 of [26], for sufficiently small e, 

y(a;) = 0(e^*^+'l^l) as 2; ^ -00 
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Letting, W{x) = e"^0c(a;)"yi(a;), 

M^(a;) = 0(e(^*+'')^+"l"^l) as a; ^ -oo and W{x) = 0{e^>'^+''^'') as x ^ +00 

From (|3.36p . /ii + a < < /i* + a, so 14^ will decay exponentially fast at ±00. Hence it is an solution to 
the eigenvalue problem AaW = XW. 

Conversely, if we have an eigenfunction, then it must satisfy 



W{x) = 0(e(^i+°)^) as a; -> +00 and W{x) = o(e('^i+"'^) as x 



-00 



Y{x) = e~°'^W{x) will then satisfy (|3.23p in a classical sense, although it may not be in L^. However it will 
satisfy the necessary decay estimates on y, constructed from Y as in (|3.24p . such that D{X;j) = 0. This 
concludes the proof of part (b). 

The proof of part (c) follows that of Lemma 2.9 from [27]. First, it is proved that if, for a given 7, A is 
a zero of order k of D(A;7), then A is an eigenvalue of Aa of algebraic multiplicity at least k. It is then 
shown that it cannot have algebraic multiplicity greater than k. We omit repeating these details. Part (d) 
is then a consequence of (c) and the calculations in Appendix [Pl that D{Q;j) — d\D{Q;j) = for all 7. 



Remark 3.10 For 5RA < with D(X;^) = 0, it is likely, but not proved, that X is not an eigenvalue of 
A. 



3.3.3 The Evans Function in the KdV ScaUng 

We now introduce D^{A; 7), the Evans function for (|3.23p under the KdV scalings introduced in Section 



^^jx, X^cA-f^ = l + ^-C/(e(x);7) 



n-l 



The eigenvalue problem is now 

= A 



1 



u 



Y 



71-1 



U] d, 1 



n- 1 



-U 



Y 



Recall that U = C/(^; 7) is the solution of ((TTS)) . and for 7 = 0, C/ = the KdV soliton. 
We can construct a dynamical system formulation of (|3.38[) . defining the vector Y as 



Y = 



\ n — 1 



1 + Y 

n — 1 ' 



L^Y 



+ 7^a(i + ;^[/)"9, {l + ^,u) 



Y 



(3.37) 



(3.38) 



(3.39) 
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which satisfies 



Y = i3.(e,A,7)Y 





(2 \ m—n — 1 / 2 \ ~" 



V Afl + ;^f/(f;7) 



(3.40) 
(3.41) 



As ^ — > oo, the matrix is 



Br (A, 7) 



1 

1 -1 
A 



which has the characteristic polynomial 

P-,(z/; A, 7) = _ j'^Aiy^ -v + K 



(3.42) 



(3.43) 



A few remarks about the scaled problem. The assumptions stated in Theorem 13.41 remain the same, 
except now the matrix under inspection is B*. with eigenvalue parameters (A, 7). Y and y are related: 

y(^) = iy2{i{x)) 

At 7 = 0, is 

d^LoY = [-dlY + (1 - U^ii; 0)Y] = AY 
the KdV eigenvalue problem, (|3.19p and 



i?.(e,A,0) 



1 \ 

-1 

A 



is the matrix for the KdV dynamical system. 
Proposition 3.11 (Scaled Evans Function) 
(a) I?^(A;7) is defined and analytic on the set A C C^, 

A = {(A, 7) I 7 e [0, 1) and A € A^} 

where, for 7 > 0, 



(3.44) 



A^ = <^ A I 5RA > 



C7^ 



^ I Acut(7) 



127 



27 



= {A\m> -7-^/- + 0(7-^) \ -^-^j +oi^-')-d- + oiY) 



16 



16 



27 



Ao and Q as defined by Theorem \3.b\ When 7 = 0, 

Ao = Axdv - C \ (-(», -v/4727 
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(b) For fixed 7 G (0, 1) and A G A^, 

(c) For A e Ao, 

i?,(A;0) =i?Kdv(A) 

Proof; The proof of these statements foUows that of Proposition 2.8 in [20] and Theorems 4.9-4. llof [28]. 
Proof of Part (a): For 7 > 0, as in Theorem 13.61 we must identify a set in in which the hypotheses of 
Theorem 13.41 are vahd. Parts (i-iii) are obvious as the soUtary wave U{^;j) decays exponentiaUy in ^ and, 
for fixed will be analytic in 7. We are left to verify part (iv). The characteristic polynomial of B"^ is 
(|3.43p . As noted in [20], the roots of are realted to those of P, (|3.35p . by 

^(A,7) = 7i/(A,7) 

So Pi, will have a unique root of minimal real part for a given A and 7 when P has such a unique root for 
A — cA'y^. Therefore, for 7 G (0, 1), if A is in the set 

C7'^ 

(iv) will be satisfied. If 7 = 0, ()3.40|) . ()3.4ip coincides with the KdV system, for which Aq = C \ 



(—00, —^4/27]. Clearly, as 7 — > 0, A-y limits to Aq. 

Proof of Part (b): From part (a), A — cAj^ £ fi^ and, by construction, 

yi{x, X,-f) ^ e^^^ as X ^ +00 
Yi{^,A,j) e"'^ as^^+c5o 

At ^ = 7a;, /ii — 7t^i, A — cA'y^, yi{x,X,j) — Yi(^,A, 7). Using the Transmission Coefficient interpre- 
tation of the Evans function, we then have 

y(a;, A,7) - i:»(A;7)e''i'^ (1, /ii, X/mf as x ^ -00 
Y(e,A,7)^I?4A;7)e'^i«(l, z^i, A^f as ^ ^ -00 

implying 

DicA-f^;-f)^D,{A;-f) 
Proof of Part (c): Trivially, when 7 = 0, this is the KdV Evans function problem exactly. 

■ 

Proposition 3.12 Let £2 G (o, ^4/27^ and M > 0. Set 

O = {5RA > -£2, |A| < M} 
Then for all 7 sufficiently small, O C A-^ and 



lim D^A; 7) = D^A; 0) = Z?Kdv(A) 
7— >o 



with uniform convergence for A G C 
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Proof: Clearly, for 7 sufficiently close to zero, O C A^. Since O is compact and Di, is analytic in both 
arguments, 

snp\D,{A;j)~ D,{A;0)\ 
Aeo 

may be made arbitrary small by taking 7 sufficiently close to zero. 



Theorem 3.13 There exists 7^ £ (0, 1) such that for all 7 G (0, 7*], a < 0^,(7), there exists e — £(7, a) > 0, 
such that: 

• The only eigenvalue of Aa with SRA > — e is X = 0, with algebraic multiplicity two. 

• The only zero of D(X; 7) with 3fJA > —s is X ~ 0, a root of order two. 

Proof: Following the Proof of Theorem 2.1 in [20] we break a half-plane into two parts, a bounded set 
about the origin and its unbounded complement. The operator estimates from Proposition 13.31 will rule out 
eigenvalues in the unbounded part, and the Evans function in the KdV scaling will control eigenvalues in 
the bounded part. 

Applying Proposition [33] (a) with ei = 1/4, there exists M > 0, such that for any c G (n,4n], a < 0^(7), 
and A G flu 

nu = jsRA > ~^acj^, \M > cMj^ 

||C(A)||i2 < 1. By Proposition 13.21 (c). such A cannot be eigenvalues of Aa. Additionally, by Theorem 
EH part (b), i:>(A;7) 7^ on this set. 
Let 62 = J and set 

O = {A : 3fiA > -£2, |A| < M} 

Define 

TO = min{|DKdv(A)| min A G dO} 
By Proposition l3.12l there exists a 7* < 1/2 such that for all 7 G [0,7^,], 

\D4A; 7) - D,{A; 0)| = {D^A; 7) - r'Kdv(A)| < m for aU A e dO 

By Theorem 13.51 the only root of I?Kdv in O is at A = 0, with multiplicity two. Applying Rouche's 
Theorem (see [16] for example), for all 7 G [0,7*], 7) and -DKdv have the same number of roots in O. 
By Proposition l3.11l (c). if 7 G (0,7*] and A G O, L'*(A;7) = D{cAj^;"f); therefore on the set fie = C7^C, 
D{X; 7) also has only two zeros. Theorem 13.61 (d^ asserts that A = is a root of multiplicity two. 

If we now set 

e(7,a) = min|iac7^ i*^^^'^^} (^-^^^ 
then by Theorem 13.61 (c), A = is the only eigenvalue of Aa with 5RA > —e. 
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(a) 



(b) 



Figure 4: In (a) we evaluate D{-; c = 4) on a portion of the strip 5RA = —1/5. In (b) we evaluate D{-; c = 8) 
on a portion of the strip 3fiA = —1. Since D{X) = D{X) we only compute 3A > 0, and then reflect; this is 
the dashed curve. Both curves wrap around the origin, marked by x, twice. In both cases, n = 3 and m = 0. 



For 7 > 7^,, one may compute the Evans function numerically to assert than a given A is not an eigenvalue. 
Moreover, the winding number of the image of the Evans function evaluated on the line 5RA = — Ai < 
equals the number of zeros in the set SRA > — Ai. Therefore, one might evaluate the Evans function on such 
a line and exam the plot, as we do in Figure 21 

These plots indicate that A = is the only zero in the closed right half-plane for the values of c, n, and 
m under consideration. Up to the acceptance of these numerics, this extends Theorem 13. 131 Note that we 
do not evaluate out to — Ai + too, but merely compute at sufficiently large values of A such that we are near 
the asymptotic value of the Evans function. It may be proven that there exists Dodl) such that 



A further numerical computation will reveal that this value is nonzero. 

In the Hamiltonian case, n + to = 0, an analytical result is possible for 7 > 7*. The linearized operator. 



Corollary 3.14 For 7 G (0,7*], dcAf[(t)c] ^ 0. 




3.4 The Evans Function outside the KdV Scahng 




28 



A, may be written as A = JcLc, 

Jc = [I- (CO)]"' Jc = -Jc (3.46) 

This structure permits an extension of the Theorem 13.131 beyond the KdV regime, given below in Theorem 
13.201 However, this is absent for general n and m. 

Remark 3.15 This section is the only place where the analyticity of the Evans function in the 7 argument 
is used. In turn, this is the only place requiring analyticity of (f>c in c from Corollarv \2.10\ (b). For the results 
in the preceding section, joint continuity of D{X;j) in its two arguments is sufficient. 

Lemma 3.16 In the Hamiltonian case, the following are equivalent for SfiA > 0, a < 0^(7); 

• X is an eigenvalue of A. 

• D{X;j) = 0. 

Proof: For 5RA > 0, we know that /i* = min {3fi/X2, SR/ia} > 0, so an eigenfunction, having submaximal 
growth at — CX3, must decay exponentially fast. This is very similar to the relation in the case of the weighted 
operator from Theorem 13.61 (bV 

For 5RA — 0, the proof relies on the JL structure of the operator A. See the proof of Theorem 3.6 in [26]. 

■ 

Lemma 3.17 If X is a nonzero purely imaginary eigenvalue, then d\D(X;j) 7^ 0. 

Proof: The proof is by contradiction. Let Y+ be the corresponding eigenfunction, AY'^ = Ay+. It may 
be proven that the subspace y — span |^^, ^^^^^-'^^^ 

{LcU, v) = for all u,v ^y. 

As 3^ n ker(_Lc) = {0}, we may apply Lemma 3.3 of [26] to conclude dim 3^ < 1, a contradiction. See Lemma 
3.3 from [27] for details. 

■ 

Lemma 3.18 (Analytically confirmed for n = 2) 

In the Hamiltonian case, assuming dcJV[(t)c\ 7^ for all c, then for all c, there are no eigenvalues with 
3?A > 0. 

Proof: By Theorem 13.131 the result holds for c < c^,. We argue by contradiction to extend it beyond c^. 
Assume for some 70 > 7^, there exists Aq, 3fiAo > 0, such that -D(Ao; 70) = 0. If 3Ao ^ 0, then D{Xq- 70) — 
and there would be two eigenvalues in the right half-plane. But by Theorem 3.1 of [26], A — J^L^ has no 
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more eigenvalues (counting multiplicity) with diX > than Lc does with 5RA < 0. As is discussed in [39], Lc 
has exactly one negative eigenvalue; therefore, this is a contradiciton, so Ao is real. 

Because the number of zeros in the right half-plane is at most one, counting multiplicity, we know 
9a-D(Ao;7o) 7^ 0. Applying the implicit function theorem, we get an analytic function, A(7), defined in a 
neighborhood of 70, such that A(7o) = Aq and _D(A(7);7) = 0. 

Let O be the maximal domain of analyticity of A(7). In a sufhciently small neighborhood of 70, 5RA(7) > 0. 
For real-valued 7 in this neighborhood, we must have, by the above argument about complex-conjugates, 
that 5 A (7) = 0. Considering the power series expansion of A (7), about 70, A (7) will be real- valued for 
real- valued 7 G C 

Let 

71 = inf {7 € (7*/2, 70) n O I i?(A(7); 7) = 0} 

For all 7 €E [71, 70)7 we must have A(7) > 0. Suppose not. Then, by continuity, for some 7, we must have 
A(7) = 0. But this would imply that A = was a root of multiplicity three, contradicting the assumption 
on J\f, which ensures it is a root of multiplicity two. 

A(7) may be analytically continued down till at least 7*/2. If not, then 71 > 7*/2 and d\D{X{'~fi); 71) 7^ 
since this root must be simple. Therefore we could apply the implicit function theorem again, and extend 
A (7) below 7i, contradicting its minimality. 

But then _D(A(7*);7i) = 0, and A(7*) > 0, contradicting Theorem 13. 131 

■ 

Remark 3.19 An analogous result may be found in Theorem [26] for gKdV, gBBM, and a Boussinesq 
equation. However, the argument there is very different because it may be proven that D{X) ^ 1 as |A| —> 00. 
This does not hold for the Evans function associated with (jl.2p . due to the appearance of a nonlinearity in 
the dispersive term. 

Theorem 3.20 (Analytically confirmed for n — 2) 

In the Hamiltonian case, assuming dcJV[4>c\ 7^ for all c, Theorem \3.13\ may be extended to all 7 G (7^,, 1), 
except for a discrete set whose only possible accumulation point is j — I. 

Proof: By Lemma [3.181 if Theorem 13.41 were false for some 7 > 7^,, it would be due to a zero appearing 
on the imaginary axis. We will prove by contradiction that the set 

E = {-f e [0, 1) I there exists (3 > such that D^{iP; 7) = 0} (3.48) 

has no accumulation points. We consider only positive f3's because if i(3 is a root then so is —if]. This follows 
the proofs of Theorem 3.6 of [27] and Theorem 2.1 of [20]. 

Assuming E has a limit point, there exists a sequence, jj G E, 7j — > 70 G as — > 00. Taking a 
subsequence if neccessary, 7^ and 70 are bounded away from 7 = 1. We will now rule out large eigenvalues, 
and then argue by contradiction to rule out small eigenvalues. 
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Applying Proposition 13.31 to this range of 7 values, there will exist M > 0, such that the corresponding 
f3j > must satisfy Pj < M. Taking a subsequence if necessary, /3j /3o, Po < M, Di,{i(3j]^j) — 0, and, by 
continuity, Di,{iPo;jo) = 0. 

Note that /3o 7^ 0. A = is always a root of order at least two. The assumption dcAf[(pc] « d^D{Q;^) = 
{c'^j^)~'^d\Dj,{Q;j) 7^ forces it to be a root of order exactly two. But (3o — would imply that it was a 
zero of at least four, a contradiction. 

Applying Lemma |3.17[ di^D-i,{i[3j]^.j) 7^ and dADi,(if3o',^o) 7^ 0- By the implicit function theorem, there 
is an analytic function Ao(7) defined in a neighborhood of 70, such that Ao(7o) — i/3oj a-nd -Di,(Ao(7); 7) = 0. 
By considering the power series expansion of Ao(7) about 70, we see, by taking 7^ sufficiently close to 70, 
that Ao(7) is purely imaginary for real 7 in its maximal domain of analyticity, O. 



For all 7 € [71,70], we must have QAo(7) > 0. If not, then by continuity would exist 7 € (71,70), for which 
Ao(7) — 0, yielding a contradiction as before. 

Suppose 7i > 0. 3Ao(7i) > because, if not, then by continuity there would exist 7 S (71, 70), for which 
^0(7) = 0, leading to a contradiction again. Therefore, we may be sure that 9a^*(Ao(7i); 71) 7^ 0. We 
may then apply the implicit function theorem, allowing us to continue Aq below 71, another contradiction. 
Therefore 71 = 0. But then 3Ao(0) > and = L'*(Ao(0);0) = L'Kdv(Ao(0)), a contradiction. 

■ 

Remark 3.21 This result is limited by our inability to analytically evaluate the functional A^[(/)c]. The 
authors were similarly stymied in [39], where the orbital stability of the solitary waves relies on proving 
dcJ^[4ic] > 0. Here, as there, one may numerically evaluate the functional and observe its monotonicity in 
the speed argument. See [39] for the case n — 2. 

This result, along with ()3.8p and Theorem 13. 13|. completes the proof of Theorem 13. II 
3.5 The Generalized Kernel 

Proposition 3.22 Let c > n, a < a*(c) and assume dcAf[4ic] 7^ 0. Define 
li = dx4'c 

6 = 9c0c 



Let 



71 =inf{7e [0,70) no I Z?.(Ao(7);7) -0} 



771 - e [/ 




dx 




-1 
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For j = 1, 2, set 

= e^y^^ and r], = e'^yfij 

Then {^1,1^2} and {r]i,ri2} are biorthogonal bases for ^CTg{Aa) andkeig^A*), {^i,r]j) — Sij . They satisfy the 
relations 

Aa^l = ^^772 = 

Aa^2 = -£.1 Allji = -772 

Proof: It is easy to verify A£i = 0, A(_2 = —£,1 and A* 7)2 = 0. 

ker(L:) =span{C"-"9,0,} 

The kernel is orthogonal to 772 because 772 is an even function while the kernel is odd. Threfore, 771 is well 
defined and A* 771 — —772. 



4 Semigroup Decay 

The following result proves the convective stability of solitary waves under the linearized fiow. As we will 
do in Section [51 this may be employed to prove full nonlinear stability. 

Proposition 4.1 Assume 7 e (0, 1), a < a*(7),, and there exists £ > such that X = is only eigenvalue of 
Aa with diX > —e. Then the initial-value problem 

Wt = AaW 

w{0) ^waeH^n keigiA*)^ 
has a unique solution w{t) — e'^"*wo £ Co{[0,oo); H^) with 

\\w{t)\\H^ <Ce~''*\\wo\\H^ (4.1) 

for some C > and b > 0. 

Remark 4.2 There exists a 6max > such that (14. ip will hold for all b e (0,&rnax)- In particular, for 
7G (0 

:7*]7 ^max 2i £; £ and "/i, as defined in Theorem \3A^ 
Proof: This is based on a result due to Priiss, [29]: 

Theorem 4.3 Let B be the infinitesimal generator of a Co semigroup on a Hilbert space Z . Let b > 0. If 
there exists M > such that 

\\{XI - B)-^\\z^z < M foraimX>~b 

then We^'Wz^z < e-'". 
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Following the approach in [20] for Proposition 3.1, we first show that Aa is the infinitesimal generator of 
a Co-semigroup on H^. Examining the Fomier symbol of A'^, A'^ is such a semigroup. As equation (jB.16|) 
shows, Aa ~ A^ is a bounded operator, so we may apply Theorem 3.1.1 of [24], that bounded perturbations 
of infinitesimal generators are also infinitesimal generators. 

Consider the Hilbert space 

Z-i/inkerg(A:)^ 
equipped with the tl^ norm and the operator 

B = Aa\z 

the restriction of Aa to Z . B inherits from Aa that it is the infinitesimal generator of a Cq semigroup on Z . 
Also note that a{B) = o-(Aa) \ {0} by Theorem III-6.17 of [14]. 

Recall from Theorem 13.11 cross(^a) is contained in left half-plane, and all points in (j(Aa) \ (Tcss(^a) are 
eigenvalues of finite multiplicity. By the assumption on the spectrum of Aa , the spectrum of B is contained 
in the open left-half plane. 

We will now prove there exists a uniform bound on resolvent of B for 5RA > —5 for some 6 > 0. For 

Aep(s), 

{XI - B)-' = {\I - Aa)-' \z 

so 

||(A/ - B)-'\\z^z < ||(A/ - Aa)-'\\m^m 
The resolvent of Aa may be written as, 

{XI-Aa)-'^{I-C{\))-'{\I^A^)-' 

C(A) defined in Proposition 13.21 Since 5Rf7oss(^a) < < 0, we must have b < oj. Then, for SRA > —b, the 
Fourier symbol of (A/ — A'^) ' is uniformly bounded, so 

|1(A/ - < M' for some A/' > 0. 

By Proposition 13. 31 for 9fJA > — & > —^acj^, there exists R' > such that 

\\{I^C{X))-'\\m^m <2 for |A| >i?'. 

Therefore, 

WiXI - B)-'\\z^z < \\{X-Aar'\\H^^H^ < 2M' for SRA > -6 and |A| > i?'. 

For |A| < R' and 5RA > — e, B has no eigenvalues. B is a closed operator, therefore {XI — B)^' is 
holomorphic on this compact set (see Theorem III-6.7 of [14]) giving the bound 

\\{XI - B)-^\\z^z < M" for some M" > 0. 
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Hence for all "SIX > — min{6, e}, 

||(A-B)-i||z^z <max{2M',M"} 

and we may apply Theorem [ 



5 Prelude to Nonlinear Stability 

There are three results about our system needed before we prove Theorem ll.il First, we establish criterion 
for when a decomposition of into a modulating solitary wave and a perturbation is possible. Then we 
derive equations for the evolution of the parameters associated with this modulating solitary wave. Finally, 
we relate the H^-norm to the _ff^-norm of the perturbation. 

5.1 Local Existence of Decomposition and Continuation Principles 

In analyzing the weighted perturbation w, we wish to treat the nonlinear terms perturbatively, with the 
leading order behavior governed by the linear operator Aa- As noted in Proposition 13. 22[ the operator has 
a two-dimensional kernel. To prevent the appearance of secular terms, the perturbation must be orthogonal 
to kerg(yl*); this reveals how the decomposition of (f> into a perturbation and a modulated solitary wave, 
(|2.28p . occurs. This follows the strategy appearing in of [25] and [20]. 

Proposition 5.1 Let cq > n, a < a^(co) and ti > 0. Given Si > 0, there exists 5q > such that for any 
(/) - 1 G Ci([0, ti]:H^ n H^) satisfying 



sup 

t<ti 



||e<+««)(,^(-,i)-0,„(--coi + 0o))||Hi <^o, for some 00 eR (5.1) 



there exists a unique mapping 1 1—> {0{t),c{t)) in ([0, ii]; M^) such that, 

sup|6l(<) - 6*01 +sup|0(t)| +sup|c(t) - Col +sup|c(i)| < t<ti (5.2) 
Tfe[0-l,0,c] = ((/.(x,t)-0,(t)(y),7y,(y)) = O /orj = l,2 and t e [0,ii] (5.3) 

where y ^ x — c{s)ds + 9{t). 

The number Sq may be chosen as a decreasing function of ti . 

Proof: The proof is via the implicit function theorem. In this context, the Banach spaces are C^([0, ti]; iJ^J), 
and ([0, ii]; M'^), the latter space equipped with the norm 

sup|6l(t)| + sup|^(t)| + sup|c(t)| + sup|c(i)| 
t<ti t<ti t<ti t<ti 

The functional is 
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and it is in its arguments, permitting the use of the imphcit hmction theorem. 
Setting 

'Uo = {(l)co{x - Cot) - 1, 0, Co) 
we see T[Uo] = 0. The Frechet derivative at Uq is 

/ /e''(^-<="*),50(a;,i),7/i(a;-cot)) - 5c(s)ds + 6e(t)\ 
\ (e'^(^-^«*)(50(x,t),772(a;-coi)) +&(t) J 

The derivative acting on the {9,c) is (q P), (Bf){t) = f{s)ds. This block operator has a bounded 
inverse on C^([0, fi]; R^) — > ([0, ii]; M'^). By the implicit function theorem, there exists (5o > 0, such that 
if 

sup|l0(- - 9o,t) - 0co(- - cot)\\Hi < 5o for some 9q eR 
t<ti 

then there exists a mapping in C^([0, ii]; M'^), 1 {9{t),c{t)), satisfying 

sup|6'(t)| +sup|6'(i)| +sup|c(t) - Col +sup|c(t)| < 5i for t < h 
r[0(- - 9o, t)) - 1, 9it),cit)] =0 for t < ii 

Defining 9(t) ~ 9{t) + 9o, and applying the change of variables x t-^ x + 9o, we have the form given in (jS.ip . 
([Q]) and ([O]) . 

Finally, because our norms are taken as the supremum over t < ti, if 6q works for ti, then it will also 
work for any t2 <ti. This allows 6q to be constructed as a decreasing function of ti. 



Proposition 5.2 Let cq > n, a < a^(co) and assume 4>—lG ^^^([0, Tmax]; n H^), Tmax > solves (fTT 

Given 5i > 0, there exist Sq > and 5'q > 0, such that for any to e [0, Tmax), */ the decomposition of <f>, 
v{y,t) = (l){x,t) - (t>c(t){.y), y = X- J^c{s)ds + 9{t), and {9,c) G ([0, io]; 1^^); satisfies 

sup||i;(i)||Hi < (5o/3 (5.5) 

t<to 

sup|c(t)-co| <(5[, (5.6) 

t<to 

T[0- l,6i,c](i) = forte[Q,to\ (5.7) 

then there is a unique extension of {9,c) in C^{[Q,to +ti,];M.'^) for some t^, > such that 

T[c^~l,9,c\{t)=0 fort<tn + t^<T„,,^ (5.8) 
sup|6'(t)-0(io)|+sup|^(t)|+sup|c(t)-c(to)|+sup|c(t)| < 6i, te [to,to+i*] (5.9) 

Proof: This follows the proof of Proposition 5.2 in [27], although we are forced to modify it as we do not 
know a priori that a solution exists for all time. 
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Given Si, let (Jq > be the value from Proposition 15.11 with ti = ^(Tmax — ^o)- Set (j){x, t) = (j){x, t + to), 
6i = — c{s)ds + 0{tQ). Let S'q be sufficiently small such that 

||0c' — 4>co\\H^nH^ < <^o/3 for all c' such that |c' — co| < S'q (5.10) 

Then, since 

|je-(-+«0 (^(0) _ + 0,)) 11^, < \\e<+^^) (0(0) - (. + 0,)) 

+ ||e<+''^)(0.,(.+0i)-</-.„(. + 0i))|Ui 

= Hto)\\Hi + II0C1 -0coll//i < ^^0 

we have 

\\e<+'^'> ~ 0c„(- - cot + 0,)) \\h^ < \\e<+'^^^ (^(i) - <^(0)) H^. + 

As (j) is continuous in time, there exists a < ti, such that 

sup||e"(-+'^i) U{t) - - cot + 0i)) \\m < So 

Therefore, satisfies the hypotheses of Proposition 15.11 We have a unique (0(t),c(t)) with (6'(0),c(0)) = 
{Q(to),c{to)). This gives us the extension, {B(t), cit)) = (0{t - to) - 6*0 + 0{to),c{t - to)) for t e [to, to + t^] 
and ([5^ wiU hold. 



5.2 Modulation Equations 

Given that the perturbation must be orthogonal to kerg(^*), the associated constraints give a pair of ODEs, 
coupled to the perturbation, giving a complete system of three equations for the three dependent variables. 

Let P denote the projection onto kerg{A*). Assuming this space is two dimensional, we use the biorthog- 
onal bases given in Proposition 13 . 22l to define projection onto this space and its complement, 

P^{m,Hi + {m,H2 (5.11) 

Q = I-P (5.12) 

The secular terms are excised from the perturbation equation, (j2.36|) . by requiring: 

Wr=AaW + Qg (5.13) 

Pg = (5.14) 

Pw{t 0) = (5.15) 

Constraint (|5.14p corresponds to 

{V3,g)=0 for J = 1,2 (5.16) 
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These two equations govern c{t) and 0{t), completing our system. 

Defining pi{y,t) = dy(l)c{t){y) - dycjicgiy) and P2(y,0 = dc(l)c(t){y) - dccl^coiy), the derived system is: 




'^0 \{ri2.,Sa[Ci),C,9]w 

However, the right hand side still has 6 dependence. Observe, 




(5.17) 



Co 



-S'a[co, c, 6]w — cm 



emdy log (^,J + ^^i?, [(<^-^i - co(Co' 1)) '] 



W 



cm 



emc^'^^H^l^ [9,(log0,J> 



Defining 



si = cm {(j)^^ dy(t)co 



K ^dy<t)c) 



Si = n[^^^ZH-l^Da [(0-^ - co(C„^ - 1)) •] 



Si^ 



^"^H-^Da [(C'-c(Ci-l)) •]} 



and 



the right-hand side of (|5.17p may be written as 

V2-, SaW 



Bit) 



e 



{rii,TaW 

{V2,TaW 




Eauation (|5.17p may be solved algebraically so that 6 only appears on the left hand side, 

1 + {m ,Pl) - {Vl, TaW) {fji , P2) 

{V2,Pl) - {m,TaW) l + {rj2,P2)) \cj 




(5.18) 

(5.19) 
(5.20) 

(5.21) 
(5.22) 

(5.23) 



(5.24) 



B{t) = I + 0{\c{t) — cqI) +0(||w||i2); for sufficiently small \c{t) — co| + B is invertible. Thus we have 

equations for and c, closing the system for (w, c, 9). 

Remark 5.3 ((5?24| . we see that, when B{t) is inverted, the right hand side of the system is continuous 
in t. Therefore, provided c{t), 9{t), and w(t) are continuous in t, c{t) and 9{t) will actually be C^. 
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5.3 Lyapunov Bound 

Using the functional, A/'[0] defined in (|2.8p . we have 

Proposition 5.4 Let cq > n, a < a^,{co), and let 0(a;,t) be a solution to p.2p in C^{[0,T]; n H^) with, 
data 

Assume the decomposition (f)(x,t) — > (v{y,t),6(t),c{t)) exists for t < T and 

\cit)-co\ + \\v{-,t)\\Hi <5i<l fort<T 



^0 

C=Co 



Then there exist constants K and K' such that 

Wv^m (1 - K'Mh^) < k(\AM\ + \\wh. + \\w\\h 



+ \c{t) - col + \c{t) - Cop + \c{t) - cop) 
AAA = AA[<^,(t)+t,]-AA[</.,J 

Proof: Taylor expanding Af about a solitary wave with perturbation z. 

+z]= AA[</),J + (<5AA[0,J, z) + i {S'M[cbco]z, z) + 0{\\z\\l,) 

From [38], the first and second variations are 

//I — (7)1-"-™ \ 

= e-i(r72,z) 



(5.25) 



(5.26) 



(5.27) 

772 and as in Proposition 13.221 

Take z{y,t) = (j)c(t){v) - <l>co{y) + v{y,t) = (j>{x,t) - (f>c„{y)- Then 

((5AA[<^eo],^) = (e-i772,2(2/,t)) = e-i (772,<Ac(t)(y) - 0^0(2/)) + 6-^ {m,viy,t)) (5.28) 

Using the continuity of c i— s- (/)c — 1, 

e-i(772,0,(t)(y)-0eo(l/)) <^|c(t)-co| (5.29) 
The term with the perturbation, v, may be bounded by 

e-i (f/2,w(y,i)) = e-i (772,^(^,0) < e-ih2|U2||u;|u. < kMl^ (5.30) 
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Now we bound the second variation. For brevity, let 



Then 

n2 



> K\\z\\]j, + - + 2 - ^;) + ($co, i;^) (5-31) 

> K^\\v\\]J^ - K2\c{t) - cop - KMt) - col - lU {<^c,e~''y,w^) 

We would like <i>c(,e^^"^ G so that the last term may be estimated by Huijl^g- Since a < 0^(7) < ^7, we 
may do this. 

Lastly, we have the remainder term TZ[(f)coTz]. Because of the a priori bound on this may be 

estimated as 

mc„z]\<K\\z\\%^<K{\c{t)-co\' + \\vrm) (5.32) 
Combining these estimates, ([QS)) . (jOO)) . ([OT|) . and ([O^ , 

- ^IklUO < if (|AAA| + |c(t) - col + m - cof + m - col'^ + IklU. + ||u;||i.) 



6 Proof of Main Results 

Before proving Theorem II. 1[ we make an a priori estimate. 
6.1 A Priori Estimates 

Proposition 6.1 Let cq > n, a < a*(co), and assume there exists e > such that 

a{Aa) n {SRA > — e} = {0} , X = is an eigenvalue of algebraic multiplicity two. 

Let T > 0. There exists E (0, 1) and Ki, > 1 such that, if the 4>{x, i) — 1 G ([0, T]; H i/^) solves 
(jl.2p and satisfies, for t < T, 

m{(f>{x,t) > ao > (6.1) 

X 

mf<j){x, 0" - a^<j)ix, tT > /3o > (6.2) 

X 

and furthermore: 

1. The decomposition 4>{x, t) 1— > {v{y, t), c(t), 9{t)) exists for t < T 
ii. Fort<T 

Co 



V\^\ + V\HmH^ + V\<t) - Col + \e{t) - 9o\ 



^ c{t) - eit) 



+ ||i;(t)||Hi <(5, (6.3) 
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iii. The data satisfies 



then forte [0,T], 



v/|c(O)-co| + |0(O) - 9o\ + ^\AM\ + ^\\w{0)\\hi <e<5. 



^e-^\\w{t)\\Hi + ^\c{t) - col + \e{t) - 0o| 

with K = K(ao, /3o, S^,) G (0, ^max)- 



Co 



c{t) - e{t) 



(6.4) 



(6.5) 



Proof : The strategy for proving this proposition is to show that the left-hand side of (|6.5p may be estimated 
in terms of their data, (|6.4p . using (|6.3p . This largely follows the proof in [20], with a few changes. The 
need to estimate \1 — cq/{c — 9)\ in terms of the data will require use of the modulation equations, (|5.24p . to 
control 0^ and to control llwH/fi, we will need to work in r-time. Thus we make the following estimates: 
Temporal Change of Variables: First, let us assume that (5* < i. Then, since cq > 7i > 1, ^cq < c{t) < 



|co. Furthermore, this initial choice of Si, ensures 



1 ^ Co dr ^ 3 

2 - c{t) - 9{t) ^dt -2 



so the change of variables r ~ T{t) is well defined. 

Time Derivatives of Modulation Parameters: Examining (|5.24p . B{t) = I + 0{\c{t) — co|) + 0(||w||i2); so 
there exists 5e > such that for Si, < 6b, B{t) will be invertible. Therefore 



L2 



\m + \c{t)\<\B{t)-'\{\\SaWU. 

permitting the estimate 

m\ + m\<K,si 

As terms of the form l/(c — ^) will appear, we will assume that 

1 



Qi\\l^)^<K^{\c-M + \\vM\\w\U 



(6.6) 



(6.7) 



S. < 



2Ki 



1/3 



so 1^1 < |, and this quotient will be well defined and bounded. 
Weighted Perturbation: In r-time, j|w(T)||//i is 



w{t) 



w{T{t = 0)) + 



(0) 



By the semigroup decay estimate of Proposition 14. H there exists K2 > and 6niax > such that 



\\w{t)\\hi <K2e-'>^\\w{T{0))\\Hi+K: 



-6(t-s)| 



e ^'\\Qg{s)\\mds 



r(0) 



for any b G (0, 6max)- Estimating Q, 



Co 



(l^l + |c|) 



+ |c-co| + 



1 - 



Co 



Co 



(6.8) 



< K3Si,\\w\\Hi 
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We have made use of (|6.6p to control 9 and c in terms of ||w||jLfi. 
Therefore, 

||«^(r)||Hi <K2e-''^\\w{T{Q))\\H^ +K2K:,5, j e-'''^^-'^\\w{s)\\Hids (6.9) 

Jr{0) 

Defining i!{s) — e^^\\w{s)\\ui , (|6.9p becomes 

Jr(a) 

for which we may apply Gronwall's inequality to get 

\\w{t)\\h^ <i^2|k(T(0))||Hie-(''-^^^^^*)(^-^(°» (6.10) 

So for S-i, small enough, b — K2K^5i, > and we induce decay in the norm of w. In particular, suppose 
that 5^<5b^ ^b/{K2K3) and let 

b' = b- K2K3S, 

We then return to i-time, 

T - t(0) - - /* c{s)ds + - {d{0) - 0{t)) > - (CO - S,) t - 2- 

CO Jo Co Co Co 

Therefore, 

\\w{t)\\Hi < K2\\wit ^ 0)\\Hie-^' (6.11) 

with K = ^b'. We now have e'^'||w||/fi estimated in terms of the data. 

Unweighted Perturbation: Applying this to Proposition [531 we have the estimate, 

\\vmm < K (/^+ V|c(t)-co| + |c(0 - col + \c{t) - cd^/^ + ^\Mf)hJ + W^mm 

< K (/^+ V|c(t) - co|(l + 5. + 5l) + ^\\w{t)\\L^{l + 5,)) (6.12) 

< K (vlAAH + V|c(i) - col + ^RWl^ 



If we had control of \J\c{t) — co|, then we would also control ||w(i)||^i in terms of the data. 
Deviation in c from co: Estimating \c{t) — co| using (j6.6p and (|6.11[) . 

\c{t) - col < |c(0) - col + / |c(s)|ds < |c(0) - col + / Ki (|c(s) - co| + Ms)\\h^) \\w{s)\\L2ds 

Jo Jo 

t rt 



< |c(0) - Col + Ki5, / \\w{s)\\Hids < |c(0) - co| + / K2\\w{tO)\\Hie~^'ds 

Jo Jo 

< \ciO)^co\+KiK2S4w{0)\\m/K 

So we now have \c(t) — co\ in terms of data, which we rewrite as 

v/|c(t) - col < (v/|c(0) - col + v/||u'(0||hi) (6.13) 

which in turn gives 

Mmn^ < K, (v^AAH + v/|c(0) - col + Vlk(O)llffi) (6.14) 
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Deviation in 6 from 9^): As with the speed parameter, 

W) - 0o| < 1^(0) - 0o| + / \e{s)\ds < \e{0) -eo\+ f {\c(s) - co| + \\v{s)\\h^) \\w{s)\\L^ds 

Jo Jo 

< \diO) ~ eo\ + Ki5, f \\w{s)\\Hids < |c(0) - col + f K2\\w{tO)\\Hie''''ds 

Jo Jo 

< \eiO) - 9o\ + KiK264w{0)\\hi/k 
This is rewritten as 



\e{t) ~ 0o\ < K, [\eio) - 9o\ + V\\wio\\H^ 

Another Estimate on the Temporal Change of Variables: 



(6.15) 



1 



Co 



< '^r"':t'" <2(ic^coi+|gi 



Then using (HilD and ((O^ 



Co 



< K 



(v/|c(0)-co| + \\w{t)\\Hi) K < (|c(0) - col + K2\\wm\H^) 



<KJV\c{0)~co\ + V\\wm\H^ 



(6.16) 



Combining (f6TT|) . ([6T3l) . (j6T4ll . ([6T5l) . and ((6T6l) . we have dUS]) with 5^ = min{i, (5^, (5b, (5^, Jg}, isT^ 
max{^2, iir4, K5,K6, Kr}. 



6.2 Main Result 

We now prove Theorem ll.il Let cq S (".,c*], c* the value corresponding to 7* from Theorem 13.11 (b), and 
a < a^(co). 

Define T to be the set of non-negative numbers, T, such that, given cq, a and vq € n H^: 



a solution exists, (/> — 1 e C([0, T), i?^ n iJ^), satisfying 

1 



inf0(x,t) > - = ao > 

a; 4 

inf (0(x, t)" - a^^ix, t)") > J inf (<^co (2:)" - a'^co (a^)") = /?o > 

• a decomposition of </) into {v{y{x,t),t),9{t),c{t)) exists for t E [0,T) 

• (ESI) holds for t e [0,T). 

Set T* ~ supT. We will first show that there exists e* > 0, such that for e < e*, if 



(6.17) 
(6.18) 



ll^'o||ffl + IK'ollffi < e 

then > 0. This will be proved using the continuous dependence upon the data. Using Proposition 16.11 
we will then prove — 00. 
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Let Si, and Ki, be as in Proposition (HUl 

The most difficult part of the proof wiU be ensuring the persistence of (|6.3p . Consider, at t = 0, the left 
hand of that equation side may be estimated with 



LHS(t = 0) < \AU\ + \\vo\\m + \\dy^co\\m\Oo~e{0)\ 



+ V|c(0)- col + 1^(0) - 9o\ + II0CO - ^cio)\\m + |1 - -^^^ 1 (6.19) 

+ ^Je-m\\vo\\Hi + ^\\dy<j>co\\Hl\Oo - Om + ^Uco - 0c(O)lki 

There exists a choice of S' and e' such that if 

|c(0) - col + |c(0)| + \9{0) - eo\ + \m\ < 5' (6.20) 
\\M\H^nHl < e' (6.21) 

then the right hand side of (|6.19p will be bounded by i(5^. Set 6i = min <5'}. From Propositions 15.11 and 
15. 2[ let Sq, S'q, be the corresponding values for 5i. 

There exists eexist e (0, 1) such that if ||iio||//in_ffi < Cexist then 

inf + Oo) + vo{x)] > 2ao 

inf [{cP,„{x + 9o) + vo{x)r - a' {^c„{x + ^o) + fo(a^))"] > 2/3o 

By Theorem Em there exists h > and a solution in C^{[0,ti],H^ D Hi), satisfying (|6T7| and ([OS]) . 
Furthermore, we will have the a priori n Hi bound that. 



sup||0(t) - iWmnHi < 2||</>co(- + ^o) + V0- MlmnHi < 2(||0eo - 1|| + 1) 
t<ti 

Set 

ei = min jeexist, e', ^e-^'o^Joj (6.22) 
and let ||wo|| < ei- As noted, the solution exists, satisfying (|6.17p and (|6.18p . until at least ti > 0. At t = 0, 



||e<+«o) (0„_0^^(.+ 0^)11^, <e^'"\\vo\\H^ <^So 

so by the continuity of (j) in time, we have 

\\e<+Oo){c^(t)~^,,{--Cot + eo)\\H^ <So for some i2 & (0,ii) 

Therefore the decomposition exists, with the Si bound on the modulation parameters, up till t2 > 0. 
Also at t = 0, using the Si bound on the parameters. 



v/|AAA| + ^\\w{0)\\hi + VHO) - col 



1 



+ lk(0)|ki < ^5. (6.23) 



c(0)-e(0) 

All of these terms are continuous in time, there exists some e (0, ^2)2, for which this remains smaller than 
(5*. Therefore, for e* < ci, ^3 G T , and > 0. 
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Continuing to Infinity: A few more constraints on are needed to continue out to t = cx). There exists 
ea/3 such that for e < eap, if 



V|c - Col + llt^llffi < e 



then 



inf [0c(y) + v{y)] > ao 



inf [{My) + vo{y)r - (My) + vo{y)r] > Po 



Let £2 > be so smah that 



K^,e2 = min |y, \/%, ^'^"'^} ^^d set 

= min{ei,e2} 



(6.24) 
(6.25) 



Now, assume HwoHifinffi As above, for this data we will have T* > 0. Assume T* < oo. For any 

T < T*, on the inteval [0, T], the solution exists with (|6.17p and (|6.18p . as does the decomposition, and (|6.3p 
holds. 
Then 

Il0(t) - MlH^nHi < max{l,eH/o _ + ||«(t)||HinHl) 



< ga((co+5*)T. + |0ol+5*) 



sup I 

\c—co I <<5* 



and this bound is uniform in T < T*. By assumption, equations (|6.17p and (|6.18p hold for t G [O-.T], 
uniformly in T < T*, which may written as 



< — < oo and 



1 



1 

<-^<oo 



(•,t)™-a20(.,t)r: 

Therefore, according to (|2.7p . (j){x, t) may be extended beyond T, by some amount t2 > 0. Hence, if ^ oo 
it must either be a failure for the decomposition to continue to exist or for (|6.3p to hold. 
Again using the Proposition l6.ll and our choice of e*, 

\\w{-,t)\\Hi < 5o/3 and |c(t) - co| < 6'^ for all t < T, uniformly in T < T*. 



Since cf) exists until at least T* + i2 , we may apply Proposition 15.21 to extend the decomposition for some 
amount < <2 also beyond T.^, . 
By assumption, 

1 



^J\c{t) - Col + \\v{t)\\Hi < K^e^ < -e^p for i < 



Again, by continuity, these remains bounded by until some time e (0,t*) beyond T», so (j6.17p and 
(|6.18p also persist beyond T*. 

We may now apply Proposition past T,. This gives 

-^K|K,t)IU.<x..<i.. 



v/||z«(-,t)||ffi + VHt) - col + \9{t) - ^ol + |1 
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ior t <T < T^. As -y/j AA/"] is time invariant and smaller than i(5^ 



V\AM\ + V\\M;mm + V\cit)-co\ + m) - 9o\ + \l- — -^^TTtI + ^ ^*'* ^ a^* 

c{t) - 0{t) b 

for t < T, uniformly in T < T,. But all of these functions are continuous for t G [0,T» + ^3]; so for some 
^4 > 0, this expression remains bounded by S^,. This contradicts < 00. So a solution exists for all time, 
satisfying (|6.17p . (|6.18p . along with a decomposition and (|6.3p . 

Since we may then apply Proposition l6.1l for all time, we will always have (|6.5p . By virtue of ii'*e» < 5^, < 
Ss, we will be able to invert the matrix B{t) for the modulation equations, (|5.24p . Therefore \c{t)\ + \d{t)\ < 
ATee"'** and 

lim c{t) = Coo 



exists, and if define 

^Um (^e{t) + (c(s) - Coo) = Ooo 

then 

U{-,t) - <l>c^{- - Coot + Ooo)\\m < K^e + (• - Coot + Ooo) - (tyc{t) (■ - /o c{s)ds + e{t)^ Whi 

<K,e + K\c(t)-Coo\ + \\dy(t)cJ\HAO{t)+ f {c{s) - Coo) ds - 9 oo\ 







Similarly 



+ Coot - Ooo,t) - ^c^llffi < K^ee~ 



6.3 Remarks 

This proof is equally applicable in the Hamiltonian case, n + m = 0, for values of c not in the discrete set, 
E, of points for which Aa has an imaginary eigenvalue. 



7 Summary and Discussion 

We have thus shown that in the space n H^, the solitary waves are asymptotically stable. This dovetails 
with an extension of global existence to data in a neighborhood of the solitary waves. In the Hamiltonian 
case, we can extend it beyond via analytic continuation, as was done in [20], and this is analytically 
verified for n — 2, with computations in [39] suggesting it is true for all n > 1. Furthermore, to the extent 
that we will accept a computation of the Evans function as proof, our result generalizes to large amplitude 
solitary waves with c > c*. 

To our knowledge, this is the first result for which asymptotic stability is established for a conservative 
PDE in the absence of a variational principle. 
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Open problems include a weakening of the assumption of exponential decay on the perturbation. This 
might be accomplished through the use of an algebraic spatial weight, which would require the perturbation 
to decay algebraically rapidly. Yet less restrictive would be to use the approach of F. Merle and his colleagues, 
[10, 17,18,21]. However, there is a tradeoff in both of these approaches; weakening assumption on the spatial 
decay rate of the perturbation, weakens what can be proved about the rate at which the perturbation decays 
in time. 

Finally we remark that the multi-dimensional case of (|l.ip is wide open. While there is no existence theory 
for the two- and three-dimensional problems, perhaps a similar approach, of working in a neighborhood of a 
solitary wave, could be applied, proving both existence and stability. 

A Properties of Solitary Waves 
A.l Analyticity 

Here we provide a proof of Corollarv l2.8l Let us restate the crucial theorem. 

Theorem A.l Corollary ^.1.6 of [9] Suppose that f is a solution of the convolution equation f — K * G{f) 
such that f L'^ n L°° and lim^^^^^ /(•^) = 0. // the Fourier transform K of the integral kernel K satisfies 
the decay condition \K{S^)\ < ^i(l + ^2|Cl™) Z"?" some constants Ai, A2 > Q and m>l, and G is infinitely 
differ entiahle function whose domain contains the range R{f) of f , having all its derivatives hounded on 
R{f) and satisfying the condition G(0) = 0, then f,G{f) G H°^' . In addition, if G is an analytic function 
on an open set U containing R{f), G is continuous up to the boundary of dU of U and 



then there exists a constant ctq > such that f and G{f) both have analytic extensions to the strip 
{z e C : < (To}. 

Let Uc = 4>c — 1- By Theorem 12.51 and Corollarv l2.61 Uc is positive, in L°° n L^, and decays exponentially 
fast at ±00. Using equations ()2.14|) and ()2.15p . the equation may be written as 



d{dU,R{f)) > \\KU2 



(A.l) 



7^Wc + d^Uc + 




d^,F2{l + T;c) 



K + 1-t)' 
2 



dT = 



(A.2) 



Define 




(A.3) 



Taking a Fourier transform of (jA.2p . the equation is 



7'"c(0-e'"c(0 + GK)(0 = 



This becomes the nonlinear convolution equation 



Uc{x) ^ K * G{uc){x) 



(A.4) 



^(0 = ^-72 



(A.5) 
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u 

R{Uc) 




\ 

dist{dU, R{u^)] 



Figure 5: A plot of a possible domain U, and the range of Uc, R{uc)- Note that the distance between these 
sets as drawn is the distance into the left hand side of the complex plane that U extends, and that any such 
ovoid will be acceptable so long as it stays to the right of 3?z — —1. 

For purposes of satisfying (|A.1[) . let us take K{x) — aK{x) and G{z) ~ a^^G{z) for a > 0, a to be 
determined. Under this trivial scaling, Uc — K 't' G{uc). 

K satisfies the decay estimate for Theorem lA.il . G{z) will have a singularity at z = —1, but is otherwise 
analytic. The range of Uc is the finite segment 

R{Uc) = [0, Minax] 

and G is infinitely differential there, with all derivatives bounded. G(0) = 0. Hence, the first part of the 
Theorem lA.il applies: Uc and G{uc) are in H°°. 

Now, consider the set U in figure [S] In this figure, the d{dU, R{uc)) is the distance U stretches into the 
left half-plane. 




Picking a so small that the norm is less than 1, we can find a U such that the distance between dU and 
R{uc) exceeds H^f II2, satisfying (|A.1[) . and proving analyticity in a strip. 

■ 

A. 2 Continuity as a Function of Speed 

Consider the functional 

T [c, u]^dlu + ¥2(1 + u\c) (A.6) 

as a mapping from H'^ x R ^ L^. The solitary wave Uc — (t>c — ^ satisies T [c, Uc] — 0. Using this fmictional, 
we prove CoroUarv 12 . 101 via the implicit function theorem. Given a particular c > n, set u — u^. 
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Let -ffgyen ^^'^ -^even t>e the subspaces of and L^, respectively, of only even functions. Define the sets 



1^ c — c — n\ ^ 

Mo= (c 2~'''^^~) ^ 

No = iu& Hl^^ : \\u - u\\h- < ^ 1> C H^,. 



Z — ieven 



Note that set Mo is bounded away from zero and all of the functions in iVo are uniformly bounded from 
below by |. Thefore is well defined on Mq x Nq and will be a mapping on this set into Z. 
Set 

T = — [c,u] = d^F2{l+u;c) 
oc 

S= — [c,u]=dl + duF2{l + u; c) 
ou 



T and S are bounded operators on M — > ig^gn and ^?even ~* -^evem respectively. 

Let / e ^even Consider the problem Su = f. As an elliptic problem, this has a solution provided 
f± ker (S^)- Note that SdxU = 0. 5 is self-adjoint, has smooth coefficients, and is in one spatial dimension 
this is the unique element of the kernel. But is an odd function, hence / is orthogonal to it and the 
equation has a solution u satisfying a bound 

\\u\\m<K\\u\\L2 

Because the coefficients in S and the right-hand side, /, are all even functions, u{x) = u{—x) also solves 
Su = f. By the uniqueness of the solution, u = u, so u is an even function. Therefore u G -ffeven ^^'^ 
map S : H^^^n ~^ ^even is onto with bounded inverse. 

The kernel of S, restricted to u e H^^^^ is trivial, so the implicit function theorem may be applied to 
conclude the existence of a function Q : Mi -ff|ven! c S Mi c Mq, such that 

j^[c,g{c)] = o 

for all c e Ml. The mapping Q is C^. For any such c, 

dlgic) + F2igic) + i;c) = o 

the solitary wave equation. Therefore Q{c) = Uc = cf),-: — 1. and the mapping c (/)c — 1 is C^(M: H^). The 
analyticity of the mapping may be proven by checking the analyticity of J- in a neighborhood of (c, u). 

To prove continuity in n H^, the proof is similar. Fixing a < i, and taking c e {n/{l — 4a^),oo), we 
let Ma = Mo n (n/(l - 40^), oo), = NoHH^, and = Z n Ll. We must check that S : ff^^^i ^ Hi ^ 
^cven ^ onto with boundcd inverse. This is accomplished using the previous result and studying 

Sa = Dl + duF2{l + u; c) on H^^^ ^ L^^,,. 
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B Perturbation Expansions 

Here we provide some explicit calculations, including those for Proposition 12.121 
Note the expansions 

r = cj)'^ + nct>':~^v + f4^,,v]v 



dT 



and 



(B.l) 



B[a, b]u = -d^ (ad^u) + hu (B.2) 



Recall do]) 

dt(j) = cdc4>c + ^y'l^c + ^v'" + vt = -i4>c + H'^^^dy {4>c + «)" 

Using the above expansions and the solitary wave equation, cdycjjc — (t>TH^^dy (</)"), this may be expanded 
into: 

Vt = cj^THl'dy [-0^9^ ^ _ _ ^^^n-lQ^ (C"5«0c) v] (B.3) 

- edyv - cdc(j)c - edy(j)c (B.4) 

- mcbT-'vHl'dy ir - €) - 4>'cH^!dy iUcl^cM^) 

+ C^f^",'S[/«[0c, v]v, fm[<i^,, v]v]H-^dy (B.5) 

(|B.3p is a linear term. (jB.4p will decay to zero as 9{t) and c{t), the modulating parameters, approach their 
asymptotic limits. (jB.Sp is purely nonlinear in v. 
We define J^i[v; (j)c\, the term nonlinear in w, as 

- mct^T-'vHl'dy ir - €) - <l>7Hl'dy iU[<Pc, v]v) 

+ ^TH^^'B[U^,,v]v,U[c^,,v]v]H:^'dy (0") (B.6) 

- cj,l^Hl'B[n^r'v,rn^T-'v]H^! {H^ - H^JH^'dy (</>") 
+ ^TH^^'B[nrc~\ni^T-'v]H^!dy - 
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The operator S is given by: 
S 



Co, c, 0] = mco {^ca^y^<^o " '^~~^^<^ ^dy^l^t^ 

+ n |c^cT/a, - co(Co' - 1)) •] - - c(ci - 1)) •] I 

Finally, the terms making up Qi from (j2.36p : 



(B.7) 



^1 - ^1 + ^2 + ^3 + ^4 + ^5 + ^6 + Q7 (B.8) 

Qi = -fmicl^c, v]wH^^'dy{r ) (B.9) 

^2 = mC-iii;i7^-i [H^ - H^J H~'dy (0") (B.IO) 

^3 = -m(f>^~^wHl^dy (n^r + /n[0c, (B.ll) 

^4 = ^^TH^'aDa ifni^c, v]w) (B.12) 

^5 = ^TH^'aBa [fnicl^c, v]w , v]w] R-'dy (B.13) 

g, = -cl^TH^'aBa Hc-'w, mC~'^] {H^ H^JH^'dy (</.") (B.14) 

Sr = €'H^'aBa Hr'w, mc^l^-'w] H^^dy (0" - 0^) (B.15) 

The difference between Aa and A'^ may be written as 

+ cmC^^-^, - cn4>TH-l^Da [(1 - C')'] + ^M^a - i)-] (B.16) 

In the space weighted by (fidx)^™, this difference is 

- ir - nH^^Da [(1 - C"') •] + nH^^Da {^.^ ~ 1) H^^^ (C'^O 

- ni/^-i (C - 1)] (C^'O - nH^'a (C - 1) (C " 1) ^^^.^a^a (C^'O 
+ niJ^^ (C - 1) ^M^a - 1) D,H-^,^D, (C-^-) 

- niJi-^Z?, (c/.^ - 1) DaH-^\Da (C^'O + cmH^} (C-'^.c/.,-) (B.17) 
+ cmi/i;! (C - 1) - cmH^lDa {^c - 1) DaH-^\ {€'-'Mc-) 

+ cnH-lD, [C (C' -!)•]+ cni/^^i (C - 1) H'^^Da [C (C' - l) '] 

- cnH^^Da id^: - 1) DaH^^]^Da [cl>T (C' - l) '] 

C Analysis of the Characteristic Polynomial 

In this section we prove that (j3.35p . (A — c/i)(l — /i^) + n/i = 0, has a unique root of minimal real part on a 
slit half plane 

{A:5RA> -Ao}\(-Ao,-f^(7)] 
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There are two ways that this coiild be false; there could either be a multiple root or two roots with the 
same real part, but differing imaginary parts. As previously noted, we will have a unique root of minimal 
real part for A in the closed right half plane, so we need only concern ourselves with < 0. 

We will identify a portion of the domain SRA < for which there are neither multiple roots nor complex 
roots with the same real part. Note this is a stricter condition than is needed, as the polynomial could have 
a double root for some A, where the third root of has a smaller real part than the multiple root. 

Note that -P(±l) never vanishes, hence ^(/i) = is equivalent to = A, where 

i?(M)=CM+-^ (C.l) 
C.l Roots of Order Greater than One 

We start with the possibility of a double or triple root, as this is very easy to rule out. If is a multiple 
root, then in addition to R{fJ.) = A, we will also have 

dR__ 1 + _ n 

which has solutions 

2c + n± VScn + 

"^ = -V Yc 

We have ignored the roots with a + sign in front, as these will correspond to positive A. Note that they are 
all real, hence A will also be real. 
The A one gets from the root 



2c + n+ V8cn + rfl 3\/3 



is decreasing in c. The other root, 
will map to A values 

A(m-) : 



' 2c + n - \/8cn + 
2c 



- 8c2 + 20cn -rfi - d,c\Jv? + 8cn - n\/rfi + Sen (C.2) 
It can be checked that for c > n > 1, i? < R (n-). Therefore, for A > — ri(c), with 

f2(c) = 8c2 + 20cn - - 8c + Sen - n v'n^ + Sen (C.3) 

P(/ti) cannot have a multiple root. 

C.2 Roots of Differing Imaginary Part 

U Hi = a + and ^2 = oi-\- 1^2 are two roots of P, then 

i?(a + z/3i) =R(a + il32) 
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After matching real and imaginary parts in this expression, the three unknowns a,/3i,/?2, must satisfy the 
two equations 



Solving the (jC.4p for /3| in terms of a and (3i, there are two families of solutions: 

0l = Pi (C.6) 



Without loss of generality, we assume /3i 7^ 0. 

Recall that A is imaginary if and only if has a purely imaginary root, hence the condition 5RA < 
ensures a 7^ 0. Then, \G.7\ implies < |a| < 1. Furthermore, if \a\ — 1, then either the roots are conjugate 
or (32 = 0. But if /32 = 0, then /i2 = 1, which we know is not a root of 

C.2.1 Complex Conjugates 

When Pi = — /32 = /?, (jC.Sp implies that A is real and 

n(l + a2+/32) 



(-l+a2)2 + 2(l + a2)^2+^4 

which has roots f3'^, 



(C.8) 



f3^^-l-a' + !^±^^±2^ (C.9) 
2c 2c ^ ^ 



We may immediately rule out the negative root for /3 . For /? > 0, a must satisfy 



n Vn2 + 24c 2 . " V?^^ + 24c ,r^,r.^ 



Using (IC3| 



X{a) = 



4a 

Since we are only concerned with A < here, a must, in addition to (jC.10[) . satisfy 



4^ ^^■^^> 

When a > 0, (|C.11|) requires 

2 c - n 
<a^ < — — 
4c 

But i (c — n) /c < 1 + n/ (2c) — \/n'^ + 24c/(2c), so complex conjugate roots with a > are not possible with 
A in the left half plane. 

When a < 0, (jClip . imphes 

1 / c — n 
"<-2V^ 
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to satisfy (jC.lip . Consider a in the interval 



r n \/'n? + Sen 1 Ic — n 



In this interval, there will not be complex conjugate roots, as it violates (|C.10p . 

A, as a function of a, is negative and increasing on this interval. For A in the image of this interval, we 
may completely rule out complex conjugate roots. The image of this interval is 



A 



1 



\/n^ + 8cn 1 jc — n 



^ 17(c), 



2c 2c 

Hence, for > A > —Q,, one may rule out both multiple roots and complex conjugates. 
C.2.2 Non-Conjugate Complex Roots 

Consider the case of complex roots with the same real part, but imaginary parts such that \f3i\ ^ \f32\- 
Squaring both sides of (|C.5|) and plugging in (|C.7D for /3|, we get a sixth order polynomial in Pf. The roots, 
as functions of a, are 



f3f = -l-2a-a^ 
= -1 + 2a - 



/3? 



(n/c)2 +4(n/c)(l - 






8(1- 


a2) 




V((n/c)2 - 16(1 


-a2))((„/c)2„8(l-a2)(2a2- 


(n/c))) 




8(1 -a2) 




(?i/c)2 +4(n/c)(l - 


-a2)-8(l-a4) 




8(1- 


a2) 




, ^{{n/cr -W{1 


- a2))((n/c)2 _ 8(1 - a2)(2a2 _ 


(n/c))) 


8(1 -a2) 



(C.12) 
(C.13) 
(C.14) 
(C.15) 

(C.16) 
(C.17) 



(|C.12p and (|C.13p force (3i to be imaginary, hence they can be ruled out. Using (|C.7p . if (3i is either 
(|C.14|) or (IC.15|) . then /?2 — ±/?i, conjugate roots. 
In the last two cases, if Pf is to be real, then 



((7i/c)2 - 16(1 - a2))((n/c)2 - 8(1 - a2)(2a2 - (n/c))) > 



(C.18) 



If we can find, that for A in the left half plane sufficiently close enough to the imaginary axis it is negative, 
we will be done. (|C.18|) is negative at a = 0, so there exists a neighborhood of the imaginary axis, such that 
complex non-conjugate roots may be ruled out. 
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The roots of the left hand side of (jC.18[) are 

1 /n\2 



= 1 



16 



©' = ^(5-7^)(3 + 7^) (C.19) 



ollril/nl, , 

These are positive for 7 e [0, 1]. It may be checked that (|C.2ip is the smahest for aU 7. Hence the root of 
(jC.18|) such that /i wiU be closest to the imaginary axis is, 



(^ = ~l^^{l~l){i + l) (C.22) 

a larger than (jC.22[) and less than zero will yield a A that does not have non-conjugate complex roots. 

Given A, if A) is to have two roots of same real part, a, but differing imaginary part, then, by trying 
any of the last four roots for (3i (both positive and negative square roots of (IC.16P and (jC.17|) ), in (jC.4p 
the real part of A and a are related by 

5RA = m f 2 - — — ) (C.23) 

V 4(1 - a2) + (n/c) J ^ ' 

Note that 

d^_^( 32(1 ~ a'f + 12(n/c) ~ 20a^{n/c) + {nlcf \ 
da'^^y (4(l-a2) + (c/n))^ ) 

and derivative has one negative root with \a\ < 1 at 



Comparing (|Q:22)) with ((a25|) . at 7 = 1 

(lc:25]l = -1< = (101221) 

and at 7 = 

(|c:25i) = (Ic21 

In addition, one may check that (|C.25p is increasing in 7 € [0, 1] while (jC.22p is decreasing on the same 
interval. Therefore 



5 n 1 /„ , n 



-V' + T^c-ieV^^c+^UcJ ^-2^(i-^)(3 + ^) 

for all 7 e [0,1]. 
For Sft/i = a in 

(-^v/(l-7)(3 + 7),o) (C.26) 

(jC.23p will be an increasing function in a, (jC.24p is positive at a = 0. On the interval (jC.26p . the mapping 
is invertible and its image is 



(^-icy^^(7^-3f/^o 



(C.27) 
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Therefore if A has real part in the interval (|C.27p . and V{ii]\) is to have non-conjugate complex roots, a 
must lie in (jC.26[) . But such an a violates (jC.181) . and we may conclude that there are no such roots. Letting 
— Ao denote the 5RA value at the left end point of (|C.27p . 

Ao = ^ (7 + 3)'/' = ^(Y^ + (^-2^) 

Hence for SRA > — Aq, non-conjugate complex roots are not possible for P{^)- Together with (jC.Sp . 



D The Zero Eigenvalue 



In Section [3.5[ A = was identified as an eigenvalue of multiplicity at least two. Using the Evans function, 
the order of this eigenvalue may be related to the slope with respect to c of the invariant functional M[(j)c] ■ 
Using the framework from Section 13.3.21 set 







(D.l) 



(D.2) 



These are solutions to the dynamical systems 

y = B{x, A = 0, 7)y and z = —zB{x, A = 0, 7) 

Here, fii = —7. Employing the notation and formulation of the Evans function of [26], by Proposition 1.6, 
parts 2 and 3, since 

y+ = 0{e~'^^) as X ^ -l-oo and z = 0{e'^^)as x —00 

y+ and z~ are scalar multiples of and rj" , respectively. and rj" , the solutions of the dynamical systems 
satisfying 



(^^e^^ — > V as X — > +00 and 77 e — s- w as x — > —00. 



with 



B°°w+ = -7V+ 

v+ = (l, -7, of 

allowing us to define the Evans function as 



= (2c7^) ^ (07^, -C7, -1) 



D{X = 0) = r]-{x, A = 0) • C+(a;, A = 0) 
From the properties of the solitary waves, discussed in Section [2T2l there exists /3 > such that 

(f>^"^da;(j)ce^'' ^ -(3 as x^ +00 
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hence 

C"^ = — ^^y'^ and 77" = ^ z~ 
— p —Zficp 

D{0) = by inspection. Using (1.22) from [26], it is trivial to compute that d\D{Q) = 0. Taking the 
derivative of this formula, at A = 0, 

VxA^C+ - / 7j-AxC+ (D.3) 

-00 J —00 

and 77^ satisfy the ODEs 

y\ = Byx + Bxy and za = -zaS - zBa 
These problems are associated with the derivatives with respect to A of equations p.23p and (|3.29p at A = 0, 

d^L^Yx = [/ - a, (<^^5, (C™-))] ^ (D-4) 

-i^a^ZA = [/ - (C^:.-)] ^ (D-5) 

through the identifcations 

yf = Lc^A + red. ir^Y) 4') = -Zx 

For A = 0, ((3?23l) . l|3?29)) . p.4p and p.Sp are related to the generalized kernels of A and A*: 

Z= r §^dx Zx^-f {LIT' [l-r-d. {€d.-)] r §t^dx 

J — OQ 'rc J —00 J — 00 (Pc 

Then, 

Vx =77^{cdA€d.Zx), -c^d.Zx, -Zx) 

Z'-fCp 

Finally, we compute 

dlDiO) = -l^dMlM (D-6) 
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